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ABSTRACT 


This study began as a review of data and correlations for line plumes. In its attempt to 
develop a consistent set of correlations for that data set, it developed an analysis to account for the 
base geometry of the fuel source. An integral model was used for the plume to address the flame 
zone (near field) and the noncombusting plume (far field). Although point source models were 
used, the diameter or rectangular source dimensions were included in improved solutions for 
flame height and entrainment rate. The analysis also explicitly includes the effect of flame radiation 
by a loss fraction, X,, and the parameter, 

(1-X,)(AH,/s) 

erly 
Although radiation data are not usually recorded, estimates allowed correlations to be computed 
with radiation considerations. The entrainment coefficients determined for Gaussian velocity and 
temperature profiles and for constant density were found to be 0.098 and 0.091 for the 
axisymmetric and line far field values, respectively. The entrainment coefficients in the near field 
were found to be 0.22 and 0.59 for the axisymmetric and line cases, respectively. Theoretically 
developed equations for near field entrainment and flame height with empirical lead coefficients 
were found to fit a wide range of data including fuel base diameters of up to 0.5 m. The same lead 
coefficient is used for the line and the axisymmetric cases as suggested by theory. The entrainment 
and the flame height correlations are directly related by the ratio of air entrained to stoichiometric air 
required, n, which is found to be 9.6 for both geometries. The entrainment rate is given in the 


form, 
m, iz Wi Zz pe 
aap = cp) [1+20(5)] 


which does not depend on the energy release rate except to determine the position z where the near 
field equation no longer applies. This result is consist with a recent correlation of Delichatsios and 
an old concept originally put forth by Thomas. Yet the near field, flame zone, entrainment rate 
given by Heskestad in the SFPE Handbook (2nd ed.) depends on z and Ze, the flame height or in 
other words, the energy release rate. The figure that follows gives the essence of this near field 
result and shows how our present analytical result follows the Delichatsios correlations which also 
agree well with available data. 


w= 


Results are also presented for rectangular fires, and the theoretical basis for the correlation 
developed unifies data over a wide range of base dimensions ranging from a square to a line. 


Comparison of Axisymmetric Entrainment Correlations 
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ABSTRACT 


A unified analysis based on an integral approach is presented for fire plumes involving 
finite axisymmetric and rectangular sources. The analysis, using Gaussian profiles, obtains the 
best fits to experimental data found in the literature. Phenomenological constants in the theory are 
found to give consistent results in that coefficients expected to be numerically similar by theory are 
found similar among the various data sets. Thorough reviews of the literature data for line and 
rectangular sources are presented and yield consensus correlations, accordingly. The effect of 
flame radiation is explicitly included by a radiation fraction which proves to be a significant 
variable, previously overlooked in experiments. Effective entrainment coefficients for the far field 
noncombusting plume are found to be 0.09 to 0.10, and for the flame region about 0.22 for the 
axisymmetric and rectangle cases as long as D/L>0.1. For D/L < 0.1, the infinite line gives 
nearly double the flame entrainment. Generalized results are presented for entrainment rate and 
flame height in terms of single algebraic equations which span a wide range of Q* or energy 
release rate values. For low Q*% fires, the effect of diameter and line width are important and 
expressed by the theory but not enough to address dependence involving Grashof and possibly 
Froude number effects near the origin. Laminar results are also addressed. 


Introduction 

Many have studied fire plumes, with most attention on axisymmetric rather than line or 
other geometric plume sources [1-15 ]. Theoretical analyses are principally couched in terms of 
assumed similarity velocity and temperature distributions with Gaussian profiles proving adequate 
for the noncombusting region (e.g. Rouse [2]). Turbulent mixing is taken into account by an 
entrainment rule based on the axial plume velocity or its momentum to account for density effects. 
Often, idealized point source solutions (e.g. [3,8]) for noncombusting plumes are used to describe 
the results for combusting plumes using an origin offset, virtual source, to adjust the data (e.g. 
[7]). Steward [10, 11] developed integral solutions to both finite line and axisymmetric fire 
plumes which lead to rather complex algebraic results which might have discouraged their use. 
Thomas et al. [16] show that the flame geometry has a significant influence on the power law 
behavior for flame height. Delichatsios [17] used these flame geometry characteristics to derive 
entrainment and flame height results for axisymmetric fires which are consistent with data, thereby 
showing the intimacy of entrainment and flame height. One of the earliest and most comprehensive 
studies of plume from small fire sources was done by Yokoi [1]. He was first to point out that the 


axial velocity, w_, the axial temperature rise, T-T.,, and position, z, are related by the constancy 
of w_/ [(T-T,)/T,gz]. We will show that this constant differs slightly in the flame and 


noncombusting regions. 

This study will demonstrate how a theory can unify experimental correlations for fire 
plumes -- dissecting them into two regions: the “near field” combustion zone and the “far field”, 
having a point source character. An integral analysis will be used for each plume geometry and a 
second order correction will be applied in the near field solution to determine better results for 
flame height and flame entrainment. All numerical coefficients of the theoretical results will be 
determined from the available ensemble of data. We shall address finite axisymmetric, line and 
rectangular fire sources of width D and (longer) length L. Turbulent analyses will be considered, 


but we will draw on the laminar results of Roper et al. [18,19] to fully describe the data in the 


literature. 

This work was motivated by the need to synthesize the literature results for line plumes into 
a presentation comparable to the level of engineering comfort currently associated with 
axisymmetric plume formulas. Yuan and Cox [14] recently reviewed the line plume literature and 
presented some new data. We examined the available data in the literature to establish overall 
correlations with our theory. In many cases we have only plotted representative data due to the 
tedious nature of extracting data points from literature graphs. 
Model Formulation 

A point (or infinite line) source model is used for both the near and far field domains. For 
these cases, the base dimensions, D and D/L become zero. We make the following assumptions in 
the spirit of deriving dimensions relationships consistent with the commonly measured plumes 
variables. 
1. Use the Bousinnesg assumption, density constant except in the buoyancy term; and all other 


properties are constant. A correction factor in terms of T/T, could be included, but we regard the 


flame temperature factor as nearly constant. 


2. Gaussian profiles are used for wiw and (Teva in both regions for consistency in this 


presentation, although we have no justification for using it in the near field. For example, 
ee r\ 
T_.—T, = exp -B{E] (1) 


where B is selected as 1 for the near field and the rectangular source amalysis. 


3. Entrainment is based on a constant entrainment, a, which is the ratio of the entrainment to axial 
velocity. Without any loss in generality for the correlations, we do not represent the entrainment 


coefficient as a momentum ratio which would introduce (p/p). 


4. The energy release rate is assumed to be uniformly distributed in the combustion zone and 


directly determined from the entrainment rate by the quantity, AH /sn where AH. is the heat of 


combustion, s is the stoichiometric air to fuel ratio, and n is the ratio of air entrained to air involved 


in combustion. The last term, n, is solely a fluid dynamic factor characteristic of the mixing in the 


plume. AH /s taken as a constant for computations as 2.91 kJ/g, corresponding to methane. 


5. Flame radiation is accounted for by a constant fraction, X_. 


6. Mass, momentum and energy flow rates of the supplied fuel are small and ignored for these 
buoyancy dominated plumes. 
7. The plume is fully turbulent. 
Combusting Plume Structure: Axial Temperature and Velocity 

The details of the integral solutions are common in the literature (e.g.[3,8,10]) and the 
specifics for our results are found in the thesis by Grove[20]. Dimensionless variables are 
introduced as follows: 


T-T Win b Z 
»W=7_i B=, and C= 5 (2) 


c 


where the characteristic plume dimension, Z., is given as 


2U5 ; 2/3 

Tire pee Se and 2) ee OS (3) 

1/2 1/2, 
p Fou ys prceliog 


for the axisymmetric and line plumes, respectively. The solutions can be expressed in power 


relationships as 


B=C,C™, W=C_¢", and @=C,CP (4) 


The theoretical solutions and the experimentally determined coefficients are given in Tables 


1 and 2. Using the data of Yokoi[ ] for the axisymmetric far field, we determined the coefficients 


using a representative value of x = 0.2. Consequently, we find a value for a of 0.098 while 


Zukoski [21] reports 0.11 for no radiation loss. B is found to be 0.913, independent of x. The 


temperature and velocity data for the infinite line source are shown in Fig. 1 and 2 for the array of 


data found in the literature [1,2,8,12,14,15,22]. Here we assign a value of X_ = 0.3 as 
representative of these data and derive the empirical coefficients accordingly. We find C, = 2.6, 


the same as Yuan and Cox [14] and most others reported by them. In the combusting zone, the 


parameter 


aad i 


occurs in the solution. Except for X_ this is a constant for all of these data. The empirical 
coefficients shown in Table 2 are nearly identical for velocity and differ somewhat for temperature. 


The C, value here depends on the accuracy of the temperature measurement. For the axisymmetric 


case, we have used the results of McCaffrey who reported AT of about 800 *C[4], but Cox and 


Chitty [23] found the flame temperature for natural gas in square burners to be 960 to 980 °C for 


corrected themocouple measurements. Using the latter, increases the axi-C, from 0.347'¥ to 


0.427, closer to the line value, 0.450¥. Indeed, for large pool fires where smoke blockage 
causes X_to become zero, the maximum temperature rise is indicated to be roughly 1000 to 1250 
°C corresponding to our coefficients: C= 0.347 to 0.427. These temperatures are consistent with 


data presented by Baum and McCaffrey [24] for D ranging from 6 to 30 m. The velocity fit is 


shown in Figure 2, and the scatter in the data and the effect of X_ should be noted. 


Finally, we note the constancy of the parameter, W/E, being approximately 1.6 and 1.4 


in the far and near fields, respectively for both geometries. Also the a values only vary from 


approximately 0.09 to 0.1 in the far fields, but differ significantly for the near fields. The near 
field results come from the respective flame height analysis which will be discussed below. 
Entrainment 


The entrainment rates from the integral solutions follow from: 


dm, din: : 
ie Otay es or m,=1p,w,b (6a) 
dr, di, gsr 1/2 

and ip meas reg ae Orsime=T pW ob: (6b) 


The far field results follow in Table 3 using the coefficients listed in Table 1. Yuan and Cox [ ] 


find for their far field line plume data, C, = 0.51 while we find 0.54. 


In the near field we seek to find a solution that now accounts for the finite fuel source. 
Thomas [9] was first to argue that flame perimeter is key to determining entrainment into the flame, 
and Zukoski [25] needed this argument to explain his entrainment data within the continuous flame 
zone for burners of D = 0.19 and 0.50 m. We account for this fuel geometry effect by expressing 
the next order solution to b as 


D (7) 


Substituting Eq. (7) into Eq. (6), along with Ve and C, from Table 2, gives the entrainment results 


in Table 4 for the near field flame region. The form of the result is 


Zz 1/2 z 2 
panne 7 ola) [els] 


for the axisymmetric case and similarly for the line case as shown in Table 4. 


The best fit for C, comes from the following flame height results, and leads to the @ values 


listed in Table 2. These entrainment coefficients are just modeling constants used to account for 
mixing and could be adjusted for density differences; however, they may be inadequate to address 


pressure contributions to entrainment near the base of the plume as described by Smith et al. [26]. 


We speculate that the distinctly larger value of & for the line plume in the flame region is due the 


narrow strip flame enhancing mixing from distinct “fire whirl” flame structures which tend to occur 
near its base. As the plume rises, the line or rectangular plume tends to approach the axisymmetric 
case and we would expect these coefficients to become similar in the far field as shown in Table 1. 


The C, for the axisymmetric case comes from the best fit of the Zukoski data [25] shown 


in Fig. 3. It should be pointed out that this is the first unified single-equation correlation of such 
entrainment data. No comparable flame entrainment data exist for the line plume so we are stymied 


on a direct data fit for this C.. However, we note from the theory, that these coefficients are nearly 


1/2 


equal (and not dependent on entrainment), i.e. 7t/ (2V 3) =m “/2. Therefore, we select the same C, 


Re aret : 1/2 
for the infinite line as in the axisymmetric case, C = 0.0565'¥ ~. Alternatively, an attempt to use 
e€ 


the lowest entrainment data point of Yuan and Cox [14], that may have been at the flame tip, 


produced a C. nearly double the value in Table 4. 


Flame Height 


The deduction of flame height comes from Assumption 4 giving: 


m 
. 


eet * y 
CAC bee rics Out | semuerenryrs 


n(1— X ) (9) 


pz) D 


and similarly for the line case. By making this connection between Eq. (9) and Table 4 we easily 


obtain Table 5. The lead coefficient, C,, for the axisymmetric case was determined from the data 


shown in Figure 4 as presented by Zukoski [6]. The Heskestad [7] correlation is also shown for 


comparison. We selected the same value of C, for the line flames based on their approximate 
equality as given by theory. The C, values for each of the geometries were found by the best fits to 


their respective flame height data and are given in Tables 4 and 5 accordingly. The fit to the line 


flame height data is shown in Figure 5 where we plot all of the (representative) data [8-11,13,14]. 


Again we see the distinct effect of D for low oe (or Z /D)*” which changes the slope of the 


correlation. An experimental value of n can be deduced, independent of the theoretical profile 
effects, from Eq. (9) which gives C A= C. [AH /nse T.)I- Using the empirical coefficients from 
Tables 4 and 5 we obtain n = 9.6 for both plumes consistent with results by Taminini [28] and 


general literature estimates which tend to range from 5 to 15. 


Some of the difficulties for the correlations to match all of the data may be due to laminar 


flow effects. Yuan and Cox [14] have assigned their data below Z /D =7as laminar. Our 


correlation for the line case gives Z /D ~ Z /DforZ /D large, and Z JD ~ (Z MD) forZ /D 


small: a profound D effect. However, this D effect is not sufficient to account for all of the data. 
We suspect that both Grashof (Gr = gD! v’) and Froude (Fr = w(0)/(gD)"" 4) numbers are 


important factors for both laminar and turbulent domains. From the work of Roper et al. [18,19], 


we can express solutions to laminar flames as 


Le Lig Ae * 1/2 
Axisymmetric: io ie Q, Gt. =. and 
Line or slot (buoyancy controlled): (10) 
° 4/3 
Ze eee i Gri4 
py a Se 


The 4/3 power result for the line case does fit the laminar data given by Yuan and Cox [14], but the 


data listed by McCaffrey [27] of Becker and Liang (shown in Fig. 4) follow a 2 power behavior, 


not 1 as indicated in the laminar solution of Roper above. More significantly, the low (on region 


needs more attention, especially when it is realized that natural fires tend to fall between om values 


of 0.05 and 5 [27]. 


Rectangular Fire Sources 
Yokoi [1] and recently Hasemi and Nishihata [12] investigated rectangular burner fire 
plumes, with the latter study addressing larger propane fires (D = 0.1-1.0m up to 150 kW). Our 


analysis is extended to these finite sources by the following specifications: 


2 2 
1. Gaussian profiles, B =1, i.e. exp -]() + i |. 


Oe W/0C = ane a constant, which we found to be nearly identical for the line and axisymmetric 


cases with a C of about 1.5, varying for the far and near fields. 


3. Geometrical effects are addressed by representing: a = L/2 + Cz, and b = D/2 + C,z as done in 


Eq. (7), but we apply this now for both the near and far field regions. 


In the far field the energy equation yields 


Q(1-X) = 4) [p,c,w(T-T,) dx dy. (11) 
00 


By applying the conditions 1-3, we can solve for the temperature: 


Qil-X) = FCoE) [14 20(5)|[1+ 20(2)(3] 


pees (oie 
mod ~~ pei a 


(12) 


Explicitly, 


saestte el bet I8 I, 


We have chosen C, = (6/5)a = 0.132 which corresponds to the axisymmetric far field & of 0.11 


for B = 1 given by Zukoski [21]. Our fit to the data of Hasemi et al. [12], is shown in Fig. 6. We 
find C, = 2.3 or 2.9 ( 1-X)*” for X_ = 0.3 corresponding to propane. This C, is related to the 


previous far field values of the infinite line and (point) axisymmetric sources, but we have not tried 
to unify these results further. The above solution does converge to the correct limiting geometric 


behavior. 


In the combusting region, we simply fit the rectangular data as ® = 0.421\¥ which 


corresponds to the limiting values of 0.35‘¥ and 0.45'¥ of Table 3. 


W can be computed from the application of condition 2 with the selection of the constant C 


as 1.5. The entrainment equation can be written as 


co co 


m,=4 {fp ,wdxdy =1p,w ab. (14) 
00 


Accordingly the entrainment expressions can be determined for each plume region. By Eq. (9) or 


(12) where ® is constant, the flame height can be determined from the expression, 


 -20(r}(5)) 


and C, = 0.199. 


Z; 
Las GACY D 


74 1/2 
Ee oa ced 
(15) 


2/3 


1-X 


uy 


with C, = 0.00590 


We have retained the same C, used in both the axisymmetric and line cases as given in Table 5 The 
C, value for the rectangular data fit is slightly higher than the axisymmetric flame region result. At 
some point, for D/L small enough (roughly < 0.1), C, must increase to the larger infinite line value 
of Table 5. A refinement might include an inclusion of a D/L effect on C, for the rectangular case. 


Conclusions 

We have formulated a unified analysis that expresses the behavior of finite axisymmetric 
and rectangular fire plumes. The analysis gives the correct power relationships and consistent 
results for the phenomenlogical constants. It appears that for line plumes where D/L < 0.1, the 
entrainment coefficient in the combusting region is greater than that for the axisymmetric case; and 
both are greater than their noncombusting plume counterparts in the far field having a value of 0.09 


to 0.10. A systematic incorporation of flame radiation loss, X_, shows that this parameter is a 
significant factor that permeates the analysis and consequently it is recommended that future 


experimental studies report X_. Low Q’ values show distinctive dependence on D and fire source 


geometric effects which suggest the influence of Gr effects in addition to the D dependence given 


by the theory presented herein. 


Nomenclature 

a half-length variable 

b half-width or radius of plume 

B dimensionless length variable, Eq. (2) 


C specific heat ( 1.01 kJ/kg(K) 


P 
C. generic constants used in equations 
D diameter or line width 


g gravity (9.81 m/s) 


Ib length of rectangular fuel source 


m mass 


n fraction of stoichiometric air entrained 
Q energy release 

T radius in polar coordinates 

S stoichiometric air to fuel mass ratio 


ii: Temperature (T,=293 K) 


V horizontal velocity 
w vertical velocity 


dimensionless velocity (Eq. (2)) 


X horizontal coordinate parallel to length 
Xx. radiation loss fraction 

y horizontal coordinate parallel to width 
Zz vertical coordinate 


Z characteristic plume length scale (Eq. (3)) 


c 


Greek Symbols 


(04 dimensionless entrainment coefficient 
B Gaussian profile constant 


AH /s heat of combustion per unit mass of air (2.91 kJ/g) 


9 density (p_=1.21 kg/m’) 


@D dimensionless temperature, Eq. (2) 
C dimensionless height ratio, Eq. (2) 
ue dimensionless parameter, Eq. (5) 


Subscripts 


e entrained air 


f flame 
m maximum or center-line value 
fe) ambient conditions 


Superscripts 


@) signifies rate of change 

‘oY single prime (per unit length) 

QO” double prime (per unit area) 

(’’’ triple prime (per unit volume) 
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Table 1: Far Field Correlations: Coefficients = Theory, Experiment . 
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Table 2: Near Field Correlations: Coefficients = Theory, Experiment 


Dimensionless Axisymmetric (x,;= 0.20 for data ) Infinite Line (x,= 0.30 for data) 
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Table 3: Far Field Entrainment Correlations: Coefficients = Theory, Experiment 


Axisymmetric (y,= 0.20 for data Infinite Line (y,= 0.30 for data 
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Table 4: Near Field Entrainment Correlations: Coefficients = Theory, Experiment 
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Figure 1: Dimensionless Temperature vs. Height for the Infinite Line 
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Figure 2: Velocity Radiation Curves for the Infinite Line 
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Figure 4: Axisymmetric Smooth Fit of Dimensionless Flame Height, Z/D, vs. Energy Release 
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Figure 5: Dimensionless Flame Height Equation for the Infinite Line 
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Figure 6: Dimensionless Temperature vs. Height for a Rectangular Source 
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Empirical correlations are developed for both buoyancy dominated 
axisymmetric fires and infinite line fires, or fire sources having one side 
substantially larger than the other. New theory has been developed for the 
combustion zone or near field of the fire plume and the smoke plume or far 
field of the fire plume. Solutions were found for temperature, velocity, and 
plume width in both predescribed zones of the plume. Also considered were 
flame height and entrainment in both parts of the plume. Representative data 
from select authors was utilized to determine empirical constants to these 
solutions. After developing empirical correlations for the two extreme cases, 
axisymmetric sources and infinite line sources, theory was modified for 


rectangular cases and applied to available temperature and flame height data. 
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width of rectangular fuel source 
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length of rectangular fuel source 
specific heat (kj/kg-K) 

generic constant used in power series 
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mass flow rate 

multiplication factor for stoichiometric oxygen 
heat transfer rate per unit area 


energy release rate per unit volume 


energy release rate (kW) 

horizontal direction in polar coordinates 
stoichiometric air to fuel 

Temperature 

horizontal velocity 

vertical velocity 

horizontal direction parallel to length 
horizontal direction parallel to width 
vertical direction 

characteristic axisymmetric height (m) 


Z** characteristic infinite line height (m) 


Greek Symbols 
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B 
Xe: 


dimensionless entrainment coefficient 
gaussian profile constant 
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Chapter 1 INTRODUCTION 
1.1 Background 

Modeling of fire and smoke movement in buildings requires accurate 
knowledge of the structure itself including the building’s materials and construction 
method in conjunction with the potential hazard, or more specifically fuel load and 
orientation. We are addressing only one piece of this complicated matrix; the fire 
source characteristics inherent to the burning of the fuel load. Here fire source 
characteristics include flame height, entrainment, temperature and velocity both in the 
fire plume and its ascending smoke plume. Also of substantial importance is flame 
radiation which can have a dramatic effect on all of the above fire source 
characteristics as will be shown later. 

The study of fire is a fairly new field. When its earliest pioneers first studied 
and tried to quantify the behavior of a fire plume it was easiest and logical to use an 
axisymmetric source. As the study of fire progressed, although still in its infantile 
stages, scientists discovered innate differences between flat 2-D fuel sources, such as 
chemical pool fires, and 3-D sources such as wood cribs or sofas. However, studies 
were not performed to clearly show the effects that fuel geometry can have on fire 
plume source characteristics. 

As a result, numerous studies, experiments, and correlations have been 
performed or respectively developed for the above fire source characteristics for 
axisymmetric fire plumes. All pertinent information has been well documented and 
agreed upon in the fire protection community. Unfortunately, two-dimensional fire 
sources, or in this case line fires, have not received the same level of detail as their 


axisymmetric counterparts. 


Indeed, upon performing my literature review I came across a plethora of 
axisymmetric source papers that overwhelmed available information on line fires. In 
addition, a quick examination of the SFPE Handbook of Fire Protection Engineering 
addresses only the axisymmetric fire source regarding the prementioned fire source 
characteristics: McCaffrey (Flame Height) (41 Heskestad (Fire Plumes) !"°), 
Delichatsios (Air Entrainment...) [6] This is not to say that a significant amount of 
work has not been invested in the study of line plumes, merely that until recently little 
effort has been placed on comparing their correlations and thus arriving at a 
consensus for using them. Within the past few years, Gupta ?°!!?°! has compared the 
theories and analyses of some of these authors along with their empirical constants 
for the infinite line. More recently, Cox et al. ! took detailed measurements of all of 
the prementioned fire source characteristics throughout the entire plume for various 
size line sources. He then compared his empirical constants to those of more 
historical works !) 161141 

The objective of the present research is to bring the understanding of line fires 
to the same level of credibility as their axisymmetric counterparts through 
quantitative equations. This paper seeks to enhance the knowledge of two 
dimensional sources through a new theory combined with empirical constants taken 
from representative data of pertinent experiments that has been assimilated onto 
single variable graphs. Similar to the studies performed by Gupta °!?*! and Cox et 
al "! validation of previous works is also an objective. However, rather than 


examining and comparing the correlation constants of previous works on an 


individual basis the validity of these separate works will be illustrated through the 


consistency of data on single variable plots. 


Definition 


Aspect ratio is defined as the fire source’s width 

divided by length, (D/L). Theoretically, line fires have 

infinitely small aspect ratios, whereas axisymmetric D ed 

fires have aspect ratios equal to one. Aspect ratios for L 

our line data range from 0.1 to 0.007 with 0.01 being 
Figure Al: Aspect Ratio 


the norm. In other words on the average, length is 100 


times greater than the width. 


Physical Examples of Line Fires 
The importance of studying line fires and thus quantifying their behavior revolves 
around the fact that they can be used as idealized representations of situations of 
much greater complexity "7! Physical examples that correlations developed from line 
fires could be used to model include: 

e Row of townhouses 

e Long sofa 

e Advancing front of forest fire 

e Balcony spill plume 


e Flame spread over flammable wall linings PI 


CHAPTER 2 OVERVIEW OF THESIS 
2.1 Introduction 
It has been shown by Cox et al "! that a fire plume can be broken into as many as four 


separate and distinct regimes: 


e Continuous 

e Intermittent 

e Transition 

e Smoke/Buoyant plume 


Where these regimes are explained by a combination of definitions from both 


McCaffrey 51) and Morton, Taylor, and Tuner rod 


(i) The continuous zone contains a persistent turbulent diffusion flame 
yielding an accelerating flow of burning gases that are strongly 
buoyant. 

(ii) In the intermittent region temperature remains quite high, though there 
is a negligible occurrence of chemical reaction, but the gases are still 
strongly buoyant in flow. 

(i) | The transition region occurs where the velocity decays to a different 
level from that in the intermittent flame region accompanied by a 
corresponding steeper temperature decay before settling to the 
expected linear temperature decay in the buoyant plume”! 

(iv) | The buoyant plume is characterized by decreasing velocity and 


temperature with height, where the temperature has fallen to a value at 


which the strongly buoyant plume reverts to a weak buoyant plume. 


We seek to simplify the problem by dividing the plume into two regimes: 


e Continuous 
e Smoke Plume 


Some authors °!!®) have performed correlations for three regimes by inserting an 
intermittent region between the continuous and smoke plume regime. However, we 
chose to only correlate two regimes for two reasons. (1) Data can be inaccurate and 
hard to interpret in the intermittent region. (2) Individual researchers define the 


boundaries of the intermittent zone differently. 


2.2 Methodology 

Solutions were developed through a simplified version of Steward’s '*! theory. 
His analysis was diluted by incorporating point source theory; thus ignoring initial 
conditions at the fuel source including temperature, velocity, and source geometry. 
As a result of this assumption, we obtained simplified solutions that were much easier 
to use than his complicated solution. When determining plume characteristics in the 
far field this assumption proved valid as all of the characteristics fit well to available 
data. The same was true when correlating temperature and velocity in the near field. 
However, the point source assumption proved suspect when correlating flame height 
and entrainment in the near field. 


[20} were the first 


Thomas, Webster and Raftery referenced by Zukoski 
researchers to recognize the importance of source geometry, especially for large 
sources, when developing correlations for both near field entrainment and flame 
height. Therefore, through an iterative process the effects of fuel diameter were 
addressed in both of these instances where they would be the most dramatic. First, as 


stated above, the equations were-developed incorporating the point source 


assumption. Source geometry effects were then added to these “point source” 


5 


equations. The new equations fit the data much better and have a similar form to 
Zukoski’s °! corresponding correlations. 

The initial effects of velocity and plume width at z = 0 proved negligible in all 
cases. Although, Steward’s equations are more physically correct than ours the 


accuracy he gained by incorporating all of the above initial conditions is negligible. 


2.3 Outline 

The thesis can be thought of as being split into three main parts, I, II, and III. 
Part I details solutions for temperature, velocity and plume width for both 
axisymmetric and infinite line sources and in both earte of the plume; continuous, 
(near field), and smoke plume, (far field). Part II focuses on solutions for flame 
height and entrainment in both parts of the plume. However, empirical correlations 
were only developed for the near field. Part III is an extension of parts I and II and 
presents correlations for rectanglular sources normalized between their axisymmetric 
and infinite line boundaries for temperature and flame height. In general, solutions 
were developed for the axisymmetric source first as a means of verifying the new 
theory because more data is available for this source than the line. In many instances 
verification was more practical for flame height and entrainment, but nonetheless for 
consistency was carried over to the first part of the paper. 

More specifically, empirical solutions for temperature, velocity, and plume 
width are developed in both regimes for the axisymmetric source using data from 
Yokoi ®! in the far field and McCaffrey "in the near field. For the line source, the 
same sets of data [1] 7131411011021 Vere used in both regimes. The solutions to these 


characteristics are simply a constant times dimensionless height to a power. 


Temperature(®) =C,¢” 
Velocity(W ) =C,¢" 
Plume Width(B) =C,¢? 


In part IL, empirical correlations are developed for both flame height and 
entrainment. Multiple sources were used for the infinite line including Hasemi ""), 
Cox !, Steward "1 and others 5131. Although much more data 1s available for 
axisymmetric sources than for line sources only Zukoski’s °! data were used to 
develop axisymmetric empirical correlations for flame height, because we had faith in 
its accuracy and it contained a wide range of burner sizes. Zukoski’s ?°! data was also 
used to develop correlations for axisymmetric near field entrainment correlations. 
Unfortunately, no data could be found in the near field for the line so a mock point 
was invented through Yuan and Cox’s recent publication !. Unlike the correlations 
found in part I, flame height and entrainment solutions are much more complicated 
and overall more difficult to correlate. However, there is a strong inter-dependence 
between infinite line and axisymmetric correlations that was not exhibited for 
temperature, velocity, and plume width correlations. This inter-dependence is 


exploited when determining empirical correlations for the infinite line. 


2.4 Author Review 

Although the papers used in creating our dimensionless graphs originate from 
the same narrow field of study, with relatively few exceptions, the motivation for 
their study was quite different in many cases. Some authors set out to study fire 
plume characteristics explicitly NICHE PI03) Others, such as Yokoi?! and Sugawa 


(11) used line fires implicitly to model or solve new fire protection problems. Yokoi 


used an infinite line heat source in a preliminary study for obtaining the temperature 
distribution of hot currents rising from burning wooden houses in the crowded post 
World War II cities of Japan. Sugawa used line fires to study the effect downward 
vertical flow, issuing from an air conditioning system for example, would have on 
detection times of a weakly buoyant plume. Another application of line fires outside 
the field of fire science was utilized by Rouse et al '!to enable meteorologists to 
evaluate the role of the basic convective process in the more complex movements of 
the atmosphere. 

The experimental apparatus, not to mention procedures, also very extensively 
from paper to paper. Most authors used liquid fuels supplied to long metal burners 
through an externally controlled feed as their heat source. However, in many cases 
this is where the similarity ends. Burners were placed in different size rooms, with 
different ambient conditions resulting in different external wind disturbances. Some 
authors such as Lee and Emmons ™! and Yokoi ©! were careful to minimize these 
external disturbances, others were not. In addition, some burners were placed 
between parallel walls, some were set in the center of a room, and yet others were 
placed several feet off of the floor. All of these factors were noted when selecting 
appropriate data. When possible data was neglected if variances from the “norm” 
were severe. Other experimental apparatus used included long wooden strips by 
Thomas "*! and a heated wire by Brodowicz and Kierkus 7), Note that although the 
data by Brodowicz and Kierkus ware included on both temperature and velocity 


graphs it was not used in developing fire correlations. 


The following chart summarizes the data used in creating our dimensionless 
graphs. Included on the chart are the papers used; type of data recorded; type and 
amount of fuel used; and width and aspect ratio of burners. A more detailed 
explanation of the experimental apparatus and procedures used to take data is 


presented in Appendix A. 
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Chapter 3 | THEORY 
3.1 Introduction 

The objective of this analysis is to provide useful analytical relationships for 
fire plumes, concentrating on fires originating from infinite line and rectilinear sources. 
Since extensive data exist for axisymmetric fire plumes, applicability of the developed 
theoretical approach will be verified for this geometry as a justification for its use in 
infinite line and rectilinear fire plumes where supporting data may be lacking. 
Turbulent plumes shall be considered although we believe that some data exhibit 
laminar fire characteristics. 

The approach taken is based on the point source model of Morton, Taylor and 
Turner 7), and the extended analysis by Steward '! for finite axisymmetric and infinite 
line fire sources. Other theoretical analyses by Yokoi ®!, Zukoski °*!, and Lee and 
Emmons "! are also similar based on a constant entrainment coefficient. We shall 
adopt the combustion assumptions of Steward, but will ignore the finite source effects 
to allow for a direct analytical result. However, we shall, in the spirit of an iterate 
solution include the initial condition of a finite source diameter or width, but ignore 
any fuel mass or momentum source effects. It will be shown that the results produced 
in this fashion can correlate much of the combustion region effects for finite sources. 
Its implementation will also show a consistency between flame length and entrainment 
rates in the combustion zone. More specifically, entraiment rates in the flame zone can 
be predicted with knowledge of flame height. This analysis is similar in some ways to 


Delichatsios!®! and Thomas et al.°°! whose derivation is based more on flame geometry 


characteristics. Due to laminar or more specifically viscous effects, the theoretical 
results may not show all of the source geometry effects. 

In many cases a detailed examination of the solutions is only presented for the 
case of the infinite line. The interested or confused reader should find this presentation 
helpful in determining the method of solution for the axisymmetric or rectangular 
source geometries. In addition there are numerous sources available which document 


such a solution 7! 4125), 


3.2 Assumptions 

Because our methodology follows the framework laid out by Steward '*!, we 
will use all but one of his assumptions. Steward neglected radiation from the flame 
because, “its inclusion would make an analytical solution impossible.” He concludes 
that, “neglecting radiation is probably the worst assumption of all, ...because high 
temperature cases, such as in a flame, will give up a considerable fraction of their 
thermal energy by this method of heat transfer.” We, however, have included 
radiation and devoted a considerable amount of time to studying its effects on our 
correlations both in and out of the flame. 

Assumptions used in developing the theory or limitations of the theory were 
adapted directly from Steward'*! and are as follows: 


(1) Turbulent flow fully developed. 

(2) Transverse forces small compared with those in the vertical direction. 
(3) Mixing in the vertical direction neglected. 

(4) No horizontal pressure variation. 


(5) Ambient fluid is of uniform density. 


Ig 


(6) The heat capacities of the ambient air, plume and fuel are equal and 
independent of temperature. The molecular weight of the air, plume, and fuel 
are equal and the plume and air are perfect gases. 

(7) Normalized density and velocity profiles are independent of height. 

(8) The rate of entrainment is proportional to the centerline value of the plume 
velocity. 

(9) In the far field, assume an ideal point source plume where all the heat is 
released at the origin. Whereas analysis of the fire plume assumes that 
entrained air mixes with fuel and burns to stoichiometric completion 


instantaneously. 


A detailed explanation of the limitations and benefits inherent to each 
assumption can also be found in the same paper. 
In our analysis, we simplified Steward’s approach by utilizing the Boussinesq 


or small density approximation, which assumes that p = p, everywhere except in the 


buoyancy term. It can be shown that p is principally a function of temperature, T. As 


oe ni is developed. 


o o 


a result, the assumption g o ipa kot E 
3.3 Conservation Equations 

The conservation equations for continuity, momentum and energy were 
developed through a conservation of mass, force and energy in a control volume. That 
is, we define a differential control volume; identify the relevant mass, energy and force 
transfer processes; and introduce the appropriate physical equations °°. Because 
changes in the vertical direction are generally much smaller than changes in the 


horizontal r or y directions, for the axisymmetric and infinite line sources respectively, 


terms containing 0/dz are ignored compared to those containing 6/6r or 6/C y. This 
is a consequence of the far field assumption that plume height is generally much 
greater than plume width. In addition, changes in the x direction are ignored for the 
infinite line because this direction is more or less diminished by the assumption that the 
length is infinitely longer than the width. (The x-axis is taken along the length of the 
source and the y-axis is taken across the width.) A similar assumption is often used in 
heat transfer problems where the temperature is assumed constant along the axis of a 
“Jong” cylinder in a heat conduction problem. Thus a 3-D problem can be simplified to 
a 2-D problem. 

As a result of integrating the equations with respect to y or r, the laminar and 
turbulent stress terms and heat conduction flux terms drop out of the problem. The 
turbulent nature of the flow is embodied only in the entrainment term of the integrated 


mass and energy equations. 


3.3.1 Line Source 


The conservation equations for an infinite line buoyant plume are as follows: 


Continuity: Ow Ww 


—+—=+=0 
ay 
Momentum: me = Wy tae st Bs EME ot 
py a if oy 
| O'(l-xz,)_ J 
Energy: sear = it pc wT T, dy 
or 


ATs <n 6Tr \smoet paste ein Gas 
pe vZewZ)=47 X,) By 


(for combusting region) 


where QO’ = m’.AH, 


with: 
7, = The fraction of combustion energy radiated away 
mj, = The rate of fuel mass supply per unit length 


AH. = The effective heat of combustion 


3.3.2 Axisymmetric Source 


The conservation equations for an axisymmetric buoyant plume are as follows: 


ne 10 Ow 
se —=0 
Continuity =p (rv)+ a 
| pov Hw) = gp [PB 1 At) 
Momentum: /o ar 2 EP. T, Fu ey 
Asian) B 
Energy: 57 [ te oC pw(T ~ T, )ar 


or 
Chee cl swe 1 O(rg’) 
c,| v—+w— |= q4"(1- X, )-———~ (for combusting region 
poe Z =a" petro gee usting region) 
where Q = m,AH, 
with 
2, = The fraction of combustion energy radiated away 


m, = The rate of fuel mass supply 
AH , = The effective heat of combustion 


3.4 Far Field Theory for the Infinite Line, z>zp 
(Temperature, Velocity, Plume Width) 

The above equations can be written in an integral form by integrating the 
continuity and momentum equations over y from 0 to ©. Where 0 corresponds to the 
centerline of the plume and 0 physically represents the end of the domain of the plume 
or its outer radius. Upon performing these integrations, the stress term drops out of 


the momentum equation. At the plume center-line, y = 0, the stress term drops out 
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A 


ow 
because the velocity gradient equals zero, oy =0. The stress term also equals zero at 


the plume boundary, y = 00, because the ambient fluid has no shear effects on the 
plume in the vertical direction as velocity, w, equals zero. Thus the equations have 


been reduced from partial differential equations to ordinary differential equations. 


Ae ad am’ 
Cont fase. p asl CE ly yh 7 i 
ontinuity 2— J powdy = —* = -2p, lim, ,.V) (1) 
' d i.e) foe] 
Momentum: ety en (ee 
orgie) (2) 
oY b= oo) 
Energy: a ny pe,| w{T - T, Hy (3) 


We assume Gaussian profiles for temperature difference, AT, and vertical velocity, w, 


in the y- direction: 


w(y,z) _ on - y? | (4) 


w,,(z) ie 
T -T, Vi (5) 
= € _ — 
L, nm ig ta) | E b ] 
The entrainment assumption follows that proposed by Morton, Taylor, and Turner coy 
lim(-V),_,. = QW, (6) 


Where « is an empirical constant. 
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Using our Gaussian profiles, (4) and (5), in conjunction with the relationship 


fem de a vn and the substitution x = y/b the integrals in the above differential 
0 


equations, (1), (2), (3), can be solved as follows: 


(a) Jody 7 jae Pay = bf we" de = yp 


Ve lay = bw? fe de = ETE 


(b) VEEL a7 


Vi+B- b 


© folr-)o=walt, ne ay a 
T,\b 


eA ay 
21+ 8 


(d) j(r- T, dy = (T, 1 Je AZ) ah + 


_ (0-7, pvx 
2/8 


Substituting the solutions of these integrals back into their respective differential 


equations yields: 


Oana oy 


PT a5 es ee" BOR = T,)b 


Next make the equations dimensionless through the following terms: 
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ra a 


(Note- The derivation of Z** follows as the equations are made dimensionless and a 
“natural length” scale is sought since none exists for the point or line source.) 


78) =— 
(1, £8) = i (ci, os 
-" % |= 
O'"(1- 7,) = p.c, Bai Vel OTB, 


or 
esr) = :_WOB 
PL» p.¢,T,4gt? 


Let the LHS = 1 = z in order to select the most useful value to 1. Set |, 
Poe, wane i 


2/3 
equal to Z** through the following relation: /, = Z ** = (2 —o_| 
Pe pl V8 


Solve the following simplified equations a by power series: 


rade = —W (8) 


sh yeaa (9) 
de va 
1- mere le 
(l- 7,)= rio) 
Using the following assignments: 
Tha (11) 
W=C,6" 
= CLs 
After solving for the constants and powers the equations can be written: 
12 
ji 2) é (12) 
Vax 
1/6 
Br 4 ay v3 
= fh———— | 1- 
5 2B (I-x,) (13) 


_|(B a 1)" B v6(]— x). d 
oa Ty 2/3 (14) 


These constants will be determined empirically along with « and f in section 4.4 using 


two experimental inputs for temperature and velocity respectively. 


3.4.1 Entrainment 


Entrainment in the smoke plume can be computed from equation (1): 


d * dm 
20,— |way = 
Dis J in 
Recall that solution of the integral, (a) on page 8, yields: 


[vd = Yew 
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Substituting: 


d am’ 
orp ne eee 
TP, Web) 


Integrating: 
Vz p,W,,b =m, (15) 
Insert dimensionless terms, (7): 
m, = Vz p,W,b 


= J/zp,W gZ” BZ” 


= Jzp,C,\gZ"C, at 


Rearrange terms to make the equation dimensionless: 


7” 1/2 
=vacc(= ' 


= VxC,C, (07 ie (16) 


ate 


ie: 


we \3/2 oy 
Where O° = (=) - 2 
z PoC plo 82 


Empirical correlations have not been developed for equation (16) because we 
found no entrainment data other than that of Cox ”! 


3/2 


3.5 Far Field Theory for the Axisymmetric Source, z>zp, 
(Temperature, Velocity, Plume Width) 


The conservation equations for an axisymmetric buoyant plume can be written in an 


integral form by multiplying each equation by r and integrating over r from 0 to 00. 
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d * dm 
Continuity: 270,— |rwdr = —* = -27p,lim,_,..(7v 17 
ontinuity TP, zs a p, lim, _,.(7v) (17) 


a 4 r ge 
Momentum: aim dr Sali Be (18) 
Energy: O(1- 7,) t 
——_*~= p.c,|rw(T -T.)dr 
mea | G7) ag 


Gaussian profiles for vertical velocity, w, and temperature difference, AT, in the r - 


ee oa (5]| (20) 
tA a 


The entrainment coefficient, a, is defined by: 


direction are as follows: 


,s0(TV) = aw,,b (22) 


—lim 
Using our Gaussian profiles, (20) and (21), in conjunction with the relationship 
[xe de = Sand the substitution x = r/b the integrals in the above differential 
0 


equations, (17), (18), (19), can be solved as follows: 


rae) re) : 1 
rwdr = w_b*| xe dx=—w_b’ 
(e) } aia me 


r7) 0 
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¥ [rw2ar = bw? [xe ae _ Wad! 
(4) oO 4 
r yaar T,-T,’ 
(g) T -T. \dr =b?(T_-T ape Teel Ok 
[H(P—1,)ar = 0°(, ~7, fx = 


oe) * 
(h) [rw(7- T, \dr = (T., —T, w,,b? [xe P"?* de te (7, -T,)w,8" 
. 0 


2(B +1) 


Substitute the solutions of these integrals back into their respective differential 


equations: 
d 
<, wab*)= 2aw ,b 
4 (y252)- 2 Ma +1, 
dz B i% 
(7, -T, nb” _ Ol1- x,) 
(6 +1) TPC 


Next make the equations dimensionless through the following terms: 


Sarr (23) 


: 2/5 
: Q 
OO 
[2 


Inserting these dimensionless terms into the above equations and simplifying 


yields: 
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7 (WB?) = 2aWB (24) 
4 7°B')= 708 (25) 
(i-y,)= gar (26) 


Solving these dimensionless equations with the same assignments, (11), from 


B-(2a)e (27) 
5 


v-(2) Ca On a (28) 


@= 2) B'?(B+1)P% a*(1- zx, ed (29) 


section 3.4.1 yields: 


267 


Due to the large amounts of axisymmetric data and relative acceptance of 
previous correlations, comparisons of our correlations in equations (27) — (29) have 


not been made against any of these data. We rely on the data of Yokoi *!. It can be 


WONG 2m : 
shown phat a = RB is constant throughout the plume. First, Cr = 9.115 and C, 


= 3.87 from Yokoi ©! B, « and C, are then determined through a similar process 


presented for the infinite line in section 4.4. 
a = 0.098 
B =0.913 
C, =0.118 


Zukoski details a better explanation of the operation in “Properties of Fire Plumes” 


pgs. (119-121) °°! However, our values for a and C; differ slightly from his because 


24 


we incorporated a 20 % radiation loss into the equations to correct for losses from his 


alcohol data, whereas he did not. 


3.5.1 Entrainment 


Entrainment in the smoke plume can be calculated from equation (17). 


d * dm 
2ap, — | rwdr = —+ 
Po J rede = — 
Solving the integral yields: 
[rar = Orne: 
7; Z 
Substituting: 
d dm 
I w b? = e 
pL (vp) = 
Elimination of derivatives gives: 
1p,w,,b> =m, (30) 


Insert dimensionless terms: 
mM, =TP,W,0 


=2p,W.j2Z*B’Z* 
= TCG Ay A Orel ey Al 
= mp,CC; J ge *"2” 


Make the equation dimensionless: 


m 
e =aC,C; -5/6 
a anes 
=a(Gye = 
= aC, ee Ome (31) 


cee 
poVgzz? 
7 5/2 
Where Or = ¢79'? =| = 
oe Fal te: os 
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See, 


Zukoski °"! found 2C,C? = 0.21; we found a slightly lower value due to the addition 


of a 20 % radiation loss. 


3.6 Far Field Summary 


The following table serves as both a summary and comparison of the theoretical 


correlations developed for the infinite line and the axisymmetric source. 


Table 2: Far Field Summary 


Dimensionless 
Variable 


B 


” gale) 


2/3 Cala 7 iia 


> Be) 


Z**/Z* 


Et hng-x)| 6 


B 


Ba 4/3 


(2 T,Jg ) 
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(1 Je te) ih 5/3 


Infinite Line 
ea) 


(22 mea 


+1)" 
ee a a’ (1 a, 4 a FE 


par T/g | ; 


se] NI 


3.7 Combusting Plume Region, 0 < z< zy 
(Temperature, Velocity, Plume Width) 


The analysis of the combusting plume region is analogous to that of the non- 
combusting plume region presented earlier. The analysis resembles that performed 
by Steward '°! except fuel mass and momentum effects will be ignored. 

It can be shown, for a constant density Boussinesq model, that the mass and 
momentum equations are identical to (1) and (2) for the infinite line and (17) and (18) 
for the axisymmetric model. However, the energy equations for both geometries 
change in this region. 

The equations for the far field assume that all of the energy is released at a 
single point or line at the fire source. Whereas the energy equations in the 
combustion zone assume that the energy is released evenly throughout this zone and 
that entrained air mixes with the fuel and burns to stoichiometric completion 


instantaneously (6) 


3.7.1 Infinite Line 


Partial differential energy equation: 


oT 


gol Bend) )-% 


oy 


a 
It can also be shown that the energy equation in partial differential format can 
be changed into an ordinary differential equation in three steps. First, expand the two 
terms on the LHS of the equation into four terms. Two of the expanded terms then 
cancel through the continuity of mass equation. Finally, integrate with respect to (y) 


yielding: a- ae 
poe, = |w(T-T,)dy = -2,)[ q"ay 
0 0 


(32) 
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As was the case with the shear force term in the momentum equation, the heat 
conduction term also drops out after integrating. At the plume boundary, y — 0, the 
temperature of the plume has decayed to ambient, thus no heat can transfer. At the 


plume centerline, y = 0, where temperature is at a maximum, the temperature gradient 


oT Ree ra : 
is zero, —-=0. Thus, heat transfer through conduction is impossible. 


The energy release rate per unit volume integral, i} q”dy , can be related to the 
0 


flux of entrained air through: 


fara 2 3{ ay = |e (33) 


ro fa 
m’, = total rate of entraind air per unit length, 
s = mass of air required toreact with a unit mass of fuel 


However, due to incomplete mixing, s does not represent the ideal 
stoichiometric ratio. The actual air to fuel ratio is represented as ns, however 


temporarily n will not be shown explicitly but will be absorbed into s. 


_ actual air entrained 
mass of air reacted 


n 


Through an overall mass balance m’ , found in the energy release rate per unit volume 


equation, (33), can be related to density, p, and vertical center line velocity w. 
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2- | pwdy = rn! +m’ (34) 


d* n’ dm’, l 
2-— =20.aw ~<a Se 
ae Ril mm 


= the total rate of fuel flow per unit length | 


(Note- the integral in both of the above equations is multiplied by 2 to account for 


mass flow through the entire width of the plume.) 


A 


Next solve fo 


(33). 


oe {=| s Jee wdy = x - pa, 


Since s is usually large s/s+1 =1. Therefore substituting into equation (33): 


= eewaw, ca pal.) w 
n 


eat = AH,, the heat of combustion per unit mass of air reacted 
K 
AH, = 13. 1 10.2338 = 3.05-¥/ ,for most hydrocarbon fuels 
So, &. air &. air 


The full energy equation (32) can now be written: 
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—[w(T -T,)dy POW 
NCP, 
d% (T-T, (Ty AH e ® sath Sach, 
aes [ T, na ne,T, aw, ee n [way (35) 
er eee ate 
ee 


d i. 2) 
Note, through the continuity equation, (1): Yaw, = VP a he 
0 


The integral on the LHS of the equation is now solved arbitrarily assuming 
Gaussian profiles for temperature and velocity. However, in the combusting region, 
B is assumed equal to one, (B=1), since there is not enough data on the temperature 


and velocity profiles in this region to justify any other value. 


Hence: 
w T -T. —_y 
pil One y 36 
ee em”, | (36) 


Substitute the Gaussian profiles: 


Jur -T.)by = W4(Ty Te” dy = Fi =i) 
he Na 


Note, we have substitued x = y/b and used the relation fe 
, 2/2 


The integration of, (a), [way was shown to equal Vey in section 3.4, (Far Field 
0 
Theory for the Infinite Line), and will not be shown here. 


From equation (35): 


i d|¥aCa-1)5) an,» vi(E,5} 
2V2 dz ie dz| 2 


Oo 
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Make the energy equation dimensionless using the same variables presented in 
section 3.4, (7): 


Ete) how ‘es = Abe at 
ave [F=| Sony = Pa,/gZ"W = eS eZ mee) 
sefieng 


Vx d 


ee ¢ 
2/2 do 


—(W®B)= Paw =¥ WB) 


The energy equation is now in dimensionless format and can be placed with the 


dimensionless continuity and momentum equations (8) and (9) respectively. 


2a 
as a W (a) 
d¢ 7) Vx 
2 Re Ap 
rad B)=~~@B (b) (37) 


d d 
ae 18) = v2 7, 02) () 


From 37 (c) integrating with WB = 0 at € = 0: 


WOB = J2YWB 
© = /2¥ 


Plume width, B, and velocity, W, are solved for through the continuity and 
momentum equations by a power series as in the far field region, section 3.4, 
presented before. The assignments for these two conservation equations remains the 


same: 


Bei (11) 
W=C,6” 


After solving for the constants and powers the equations can be written: 
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fa fi (38) 
dad ( k 


39 
W=J¥C!” ( ) 
© = /2¥ (40) 


Empirically derived coefficients are addressed in section 4.5. The above are 
theoretical. The constant entrainment coefficient, a, is not justified, but we will see 


how well it does in fitting the above data. 


3.7.2 Axisymmetric Source 


Partial differential energy equation: 


OF ERO LS 1 O(rq’) 
pe(eZawZ) =q (i- heres 


As was the case with the infinite line, section 3.7.1, the above partial 


differential energy equation can be converted into an ordinary differential equation. 


d [2 co 
Pek p a J rw(T -T, )dr =-z,)| q"rdr (41) 


The energy release rate per unit volume integral, [ 4" rar , can be related to the flux of 
0 


entrained air as: 


farrar = {Mae (42) 
. 2a\ Sj az 
m, = total rate of entrained air 


Where s represents only the stoichiometric air to fuel ratio. As before, the 


actual air to fuel ratio is modeled through sn. Also recall that n is absorbed into s. 
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Through an overall mass balance, assuming m_,is negligible compared to m, , 
it can be shown that the mass rate of entrainment at z is equal to the mass flow rate in 


the plume at z. 


Note- m, could have been shown to be negligible through the zt = 
KY 


argument presented above. 


dm d * 
£ =27 — wrar = 279 aw_b 
= rae DOW 


Recall that < [rede = aw,_b from the continuity of mass equation, (17), found in 
0 


section 3.5. 


Next, substitute a) into equation (42): 


\ q”rdr = |= Jpnpa,b = p,aw,,bAH , 
A 2m\ Ss 


The full energy equation, (41), can now be written: 


2 [ro(7 -1,)dr _(- 7, JAH, paw b 
0 C Po 
2 trw THT, | WAZ DAH. ow b= Vow b= eS [rwdr (43) 
dz 4 T, cpl, 228 
7 AH E 


where ¥ = 
Cl, 


Recalling that # = 1in the combustion zone, the Gaussian profiles for temperature 


and velocity are as follows: 


wT Bool.) 4) 


Substitute the Gaussian profiles: 
ro) © pee (7 ™ rp 
[ro(T - 1, tr = Bw, (Z,, — T,)[ xe be = ia 
0 0 


Note- we have substituted x = r/b. 


2 


° b° : : 
The integration of wrdr was shown to equal “n in section 3.5, (Far Field 
0 
Theory for the Axisymmetric Source), and will not be shown here. 
From equation (43): 


as 2 2 
SA ee 09, Fey epi oe eae 
4 az T az 2 


oO 


Make the energy equation dimensionless using the same variables presented in 


section 3.5, equation (23): 


z* He" 


= antl Slee. a2 
=| (wos?) = Ya.) gZ° WBZ =5 VEZ" 5 qe 
oe 


d 2 d 2 
awe )=4¥aWB = 24 7 WB ) 


The energy equation is now in dimensionless format and can be placed with the 


dimensionless continuity and momentum equations (24) and (25) respectively. 


a ) = 2aWB (a) 
26 
—(W? B’*)= 
rAd )= Be (b) (45) 
d 
— (W®B’) = 2¥ — 
rai )= ae ) (c) 


Temperature,®, is solved for by integrating and canceling like terms as in section 
Sasels 
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WB? = 2YWB? | 
D=27Y 
The process for finding plume width, B, and velocity, W, was outlined in section 
3.7.1 and remains the same for the axisymmetric source. The solutions for the 


equations are as follows: 


patat (46) 
2) 

W= para (47) 

ies ids (48) 


For the same reasons presented in section 3.5, empirical correlations for 
temperature and velocity are given at this time. 
The temperature equation, (48), in empirical format: 
®=C,'¥ 
McCaffrey °*! burned methane on 0.3 m square porous burners and found that 
temperature in the flame zone was a steady 800 °C. Using this value for temperature 


and solving for ‘¥, Crr can be found easily. 


_ AT _ 800K _ 


7B 
i lg 92 Kags 

fp ETE (1~.20)2,910 _ 

On, 1.01(293) 
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Note that a value of 20 % was used for the radiation fraction. McCaffrey |"*! 


also performed extensive radiation experiments on these burners and found that x; 
varies from 10 % to 24%. These experiments are discussed in greater detail in sec. 


4.6. 


“oh 


1, =a = 0.347 
TE Ci 
® = 0.347 


® = 0.434(1- x, Y 
The empirical correlation for velocity is also solved from the same data set. 
The velocity equation in empirical format: 
= 1/2 
W= C. a4 


Choosing one representative point and changing the data presented by McCaffrey "*! 


into our format yields: 


y=1.707/ 


z=0.0726m 
Z* =0.175m 
da 
=— =0.414 
$ Z 
Next solve for W and WP: 
Wa EG 


PA ADS ae ea 


Y = 7.87 (from before) 


Thus: 
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W 1.30 


CS eee = ote = 0.720 
YE 780414 
W= 071206" 


W = 0.805(1- ,)'7C"”” 
3.8 Summary of Near Field Results 


Table 3: Near Field Summary 


Dimensionless Axisymmetric Infinite Line 
Variable 
4a 
B 4 —_—> 
rd (e k 

W ea wy Ge /2 Ve o “a 

3 
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3.9 Entrainment in the Combustion Zone, 0 < z < Z f 


The same basic equations for entrainment in the far field, found from the 
continuity equations, can be used in our analysis of entrainment in the near field or 
combustion zone. The main difference between these two analyses is the inclusion of 
geometry effects in the combustion zone. Recall from section 2.3, (Methodology), 
that both near field entrainment and flame height are developed through an iterative 
process where “point source” solutions are modified to incorporate geometry effects. 

In the far field, the role of geometry is secondary to the energy release rate 
which dominates all processes. However, close to the fire source and in the flame 
both the source geometry and energy release rate play an important role and need to 


be considered. 
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The inclusion or exclusion of geometry effects is achieved through the plume 
width term, b, found in the generic entrainment equation. In the far field, the solution 
of plume width in the form, C,¢”, is substituted into the mass entrainment equation. 
This solution only takes into account the change in plume width with height. There 
are no diameter or width effects inherent to this solution. Thus none are included in 
the far field entrainment correlation. However, source effects are inserted into plume 
width, b, in the near field through a term which includes the change in plume width 


with height, used in the far field solution, and the sources initial width or diameter; 


b=6b,+C,¢". Where 5b, is either D/2 (axisymmetric radius) or D/2 (line half width). 


3.9.1 Line Source 


Basic entrainment equation: 
m, =/rp,w,b (15) 


Terms to be substituted into equation 15: 


Where = = half width 


Remember, wm and C; are both taken from the combusting plume region, 


section 3.7.1. Substitution of the above terms into the entrainment equation yields: 


i =Apet V¥G 242" 6 


2, Se siALee ) G (=) Gy 
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Where ‘'¥ = (nx, JAF. 
c,1,ns 


In order to develop useful analytical correlations, the above equation must be put into 
an empirical format where the theoretical constants can be replaced by numerical 


constants fitted from available data. 


Empirical form: 


mm’ zV7(1 da(z 2 
— Fa = (=) (+ 6(<))2 
p,\gDD Pyke 3Vx\ DJ) 2 


rin, “ 3 BEZ a 
p, = #4) f+c(=)} (50) 


Vz 8a 


Where C, = —— — 


ee) sx 


Note- C; and C2 will be determined in section 4.11. 
3.9.2 Axisymmetric Source 


Basic entrainment equation, 3-(30): 
m, =710,W,,0° 


W = {58 vez" 


b=b, em Z ae 
2 5 


Terms to be substituted: 


Where : = radius 


Substitution of the above terms into the entrainment equation yields: 
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p.ygDD? ti 3 +$4(3)} 
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Where C, = mare CL = ee 


Empirical form: 


Note- C; and C2 will be determined in section 4.9. 
3.10 Flame Height 

Multiplying the equation for entrainment in the flame zone by a parameter 
containing the heat of combustion of air, AH, , results in a simple transition from 


entrainment to flame height. 


3.10.1 Line Source 
A useful relationship follows that energy release rate be related to the amount 
of air entrained into the flame zone times the heat of combustion of a certain amount 


of that entrained air. 


O = m'(Z,)-—= (53) 
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Using the above relation an equation for flame height can easily be deduced 


from the entrainment equations: 


es ist Spa's y Lr tn | eee 2; AH . 
P.VgDDec,T, 2 D 3Vxz\ D }J}c,T.ns 
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mies) sea aS ee 
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As before with the entrainment equations, the flame height equation is put into an 


empirical form. 


Empirical form: 


eg) 
1a D 


eWiy2 
Ties 


'(1- z,) 


Note- C and C2 will be determined in section 4.10. 


3.10.2 Axisymmetric Source 
Equality of the heat of combustion of air times entrainment and energy release rate: 


: (56) 
OQ=m, (Z,). ate 


ns 


Using the above relation an equation for flame height can easily be deduced from the 


entrainment equations: 
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Empirical form: 
Y, 1/2 Zz 2 
a5) [eolB)) ow 


p32 


'-z,) 


Note- C and C2 will be determined in section 4.8. 


DIN 


Where C = C 
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Chapter 4 EMPIRICAL CONSTANTS 
4.1 Introduction 

In this chapter, empirical constants are fitted to the theory developed in 
Chapter 3. Generally, a discussion of the data and the graphs that these empirical 
constants are taken from will proceed the explanation of how the constants were 
mathematically obtained. 

It is important to note that the data presented on all of our graphs and used to 
develop our correlations has not been weighted or selected under any mathematical 
pretenses whatsoever. Selection of data points was made in an effort to capture the 
range of data and was otherwise taken on a purely random basis. In many cases only 
limited data points were selected from a source graph so as not to overshadow data 


from other authors. 


4.2 Infinite Line Temperature Data 

Empirical correlations for temperature both in the smoke plume and the 
combustion zone are taken from Figure 1. The solid lines represent our correlations 
of the data with the flat line corresponding to the near field or combustion zone and 
the following line, having a slope of “1, representing the far field or smoke plume. 
The dashed lines depict Cox’s "! definition of the transition and intermittent zones. 
Remember that our theory models a simplified plume containing only two regimes. 
However, an actual plume can be more accurately broken into four parts, as were 
described earlier, with the addition of the transition and intermittent zones. 

To get a sense of the range of data plotted on the temperature graph lets take a 
look at the two extremes; the bottom of the flame zone and top of the smoke plume. 
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Dimensionless Temperature vs. Height for the Infinite Line 


Figure 1 
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At the bottom of the flame zone z/Z** = 0.20. In this area measurements were taken 
so Close to the bumer that they actually dip into the fuel rich zone slightly above the 
burner’s surface where combustion is just beginning to occur. Thus, these 
temperature readings are slightly lower than the other readings taken throughout the 
combustion zone. 

For large values, z/Z** = 100, energy release rates are extremely small 
compared to the heights where they are measured. Two sets of data in this area give a 
physical illustration of the large ratios between height, [m], and energy release rate, 
[m], corresponding to the very top of the smoke plume. Brodowicz and Kierkus |"?! 
measured temperatures above a heated wire with a very small energy release rate, 
9.75 x 10° kW. Thus, even measurements close to this source will yield a large ratio 
between height and energy. In the other data set, Rouse, Yih, and Humphreys were 
attempting to model atmospheric disturbances. They took measurements a couple of 
meters above recessed gas bumers emitting weak blue flames. 

Note that when developing the far field temperature correlation data points 
corresponding to large values of z/Z** were given more weight than those low in the 
smoke plume to avoid complications inherent to the intermittent zone. The same 
methodology was used when choosing the location of the far field velocity line. 
Thus, both of these far field correlations may provide erroneous answers low in the 
smoke plume. However, this area is generally not as important to the fire engineer as 
temperatures high in the developed smoke plume where interactions with detectors 


and sprinklers occur. Using these far field correlations for modelling purposes or 
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prediction of device activation is merited because they will yield conservative values 


of both temperature and the convective heat transfer coefficient, h. 


4.3 Infinite Line Velocity Data 

Empirical correlations for velocity both in the smoke plume and the 
combustion zone are taken from Figure 2. Once again, the solid lines represent our 
correlations of the data. The flat line corresponds to velocity in the far field or smoke 
plume. The preceding line, having a slope of %, represents velocity in the near field 
or combustion zone. As in Figure 1, the dashed lines represent Cox’s depiction of the 
intermittent and transition regions. 

Although the velocity data falls together in the same general area it does not 
plot as tightly together as its corresponding temperature data. This could be expected 
since it is a lot harder to take accurate velocity measurements than temperature 


measurements, therby resulting in a lack of precision among different data sets. 


4.4 Infinite Line Far Field Correlations 

Cr and Cy are determined from Figures 1 and 2 respectively by simply 
applying a best fit line through the data while maintaining the predetermined powers 
of dimensionless height from Chapter 3. These constants are then used to determine 
a, B, and C). Hence, C7 = 2.6 and Cy = 2.0 from the far field “fits” on Figures 1 and 
Be, 

Looking at either of the graphs, the constants depicted as Cr and Cy may not 
be apparent as they are expressed in a generic form to incorporate radiation loss. To 


get the values of Cr and Cy listed above a radiation fraction of 30 % has been 
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incorporated. For Cr, the radiation fraction was set at 30 % to maintain consistency 
with the propane fuel used by Hasemi |"), For Cy, 30 % was chosen for the radiation 
fraction because it did the best job of fitting available data. A detailed explanation of 
the reasons for setting x; = 30 % is presented later in the thesis in sections 4.6.1 and 
4.6.2. 

The above values for temperature and velocity agree well with the new data 
published by Cox "! and that of more historic works including Yokoi !, Rouse et al 


(41 Lee and Emmons ”!, and Hasemi and Nashimeti!"), 


2 


2 
B is first solved from the relation aa C = 2 and found to equal 0.845. With 
fC, YB 


knowledge of B, « is easily solved through the definition of Cy. 
a = 0.091 
C, is found easily as it is proportional to a. 


C; = 0.103 
4.5 Infinite Line Near Field Correlations 


4.5.1 Temperature 
First put the theoretical equation, (40), into an empirical form. 
O=C, 
Next choose a value for ® in the combustion zone that represents the majority 


of data points and solve for 


@=3.1 
ere Cael (i- (1-.30)2,910 910 _ aes, 
on 1.01-293 
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Here x, was also assumed equal to 30% to match Hasemi’s propane data |"). 
Solve for Ctr 


® 311 


Cr 


= — = 450 
sow 


6.89 _ 
Therefore: 


® = 0.450¥ 


4.5.2 Velocity 

As was the case in the previous section, put the theoretical equation, (39), into 
an empirical form. 
(ea OC Tae 

In order to solve for the empirical constant an “average” data point is selected 
and 'Y is solved for with a different radiation fraction, x,, of 15%. Here the radiation 
fraction was set at 15 % to match Cox’s methane data 
W=2.00; C=0.748 


ric) 29 ppl id) A) (Oh oh 


ef, 1.01-293 


Solve for Cy¢ 


1/2 
- _ VEE"? _ f2.89(0.748) _ tet 


a W 2.00 


Therefore: 


W = 0.800082"? 
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4.6 Radiation Effects 


The effects of radiation are considerable on our empirical correlations as was 
illustrated above. Therefore, choosing the correct x, is of the utmost importance. 
Unfortunately, x, was either not reported or not measured by any of the authors. An 
example depicting the severity of this oversight can be found on Figure 1. One would 
expect Yokoi’s alcohol data to have higher temperatures than Hasemi’s propane data. 
However, this is not the case. Had the authors reported radiation losses these values 
could have been used to correct the data. 

In most cases, the types of fuel burned were given. Although this information 
is useful it does not convey the whole picture. Radiation loss, x,, varies with both 
source size and configuration, fuel flow rate, combustion inefficiencies, and most 
importantly how the fuel is burned, or how the fuel and oxidant are mixed and ignited 
[18] 

Because x,;’S were not given we were forced to pick values based upon sound 
reasoning. Sections 4.6.1 and 4.6.2 discuss both the reasoning behind our choices and 


the effect of radiation on our temperature and velocity data respectively. 


4.6.1 Temperature 

Recall that a value of 30 % was selected for the radiation fraction, x,, in both 
the far and near fields, and was used primarily for determining a and Cyin the far 
field and Crrin the near field. The radiation fraction was set at 30 % to maintain 
consistency with the propane fuel used by Hasemi '"!. The value was determined 


from Tewarson’s chapter in the SFPE Handbook”! as follows: 
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Ey 5 1D ag Ie 


We — Frat x} 00 = 1% 100 = 28. 6% = 30% 


ead Ee Oe 43.7 

In addition, a value of 30 % gave the most reasonable radiation curves found 
on Figure 3. 

Figure 3 illustrates the effects of radiation ranging from no radiation loss to a 
hefty 60 % radiation loss. Increments of 15 % are spaced out evenly between these 
two extremes, where 30 % was taken as our baseline radiation value. As can be seen 
on Figure 3, no radiation loss will result in higher temperatures than an entity which 
radiates a fraction of its energy away. In the far field, zero percent energy loss results 
in temperatures that are at least twice as high as their counterparts loosing 60 % of 
their energy to radaiation loss and as much as 2.6 times greater in the flame zone. 

The strong effect of radiation on temperature or vice-versa was alluded to 
above. The hotter the temperature the stronger the effect of radiation. Thus, radiation 
loss in the flame zone is far more significant than losses in the smoke plume. The 
larger gap between the radiation curves in the flame zone than those in the smoke 
plume are a testament to this phenomenon. 

All of the data points within the flame zone and the smoke plume, with the 
exception of those points falling in the designated intermittent zone, are captured by 


our theory incorporating a range of radiation losses. 


4.6.2 Velocity 
Both the far and near field correlations were originally intended to be fitted to 


the natural gas data of Cox ”! and to other points belonging to low radiation loss 
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Figure 3: Temperature Radiation Curves for the Infinite Line 
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fuels. Because this LNG was composed of 94 % methane a x, suitable for methane 
was chosen. 

McCaffrey !'*! burned methane on 0.3 m square porous burners and found that 
radiation loss varies from 10 % for low energy fires to an asymptotic value of 24 % 
for high energy fires. From Figure 1 of his report we chose a value closer to 10 %. 
First of all, McCaffrey’s burners were much larger than any of those used in this 
paper. Because radiation loss increases with increasing source size!| we were 
inclined to pick a small x,to match our small burners. In order to get a more refined 
value we once again turned to Tewarson’s chapter in the SFPE Handbook 7) 
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Figures 4 and 5 illustrate the effects of radiation loss on vertical velocity. As 
was the case with the temperature plot, Figure 2, a range of radiation losses from zero 
to 60 percent was plotted with 15 percent increments between these two extremes on 
both graphs. 

Figure 5 was plotted in conjunction with Figure 4 to show the dramatic effects 
associated with choosing the baseline radiation fraction used in developing our 
correlations. Assuming an initial radiation fraction of 30 %, for both the far and near 
fields, instead of 15 % shifts all of the radiation curves upwards and captures many of 
the data points missed in the far field when assuming an initial radiation fraction of 
15%. Asa result, an initial radiation fraction of 30 % rather than 15 % was used to 
fit far field data. McCaffrey ''*! offers an explanation for this unexpected behavior: 
“Short bluish flames at low flow rates will obviously radiate proportionally less than 


the larger, yellower sooty flames of higher flow rates.” Since Cox’s "! flow rates 
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ranged from 2 — 110 kW and Sugawa (5) used propane it is posible that some of the 
data actually do approach higher radiation losses. Therefore, x, = 15 % was used in 
the flame zone and x, = 30 % was used in the far field. A short explanation of how 
radiation loss affects velocity through temperature is presented below. 

Theoretically, lower radiation loses yield higher velocities through a direct 
relation between temperature and velocity. Less radiation loss results in higher 
temperatures, Figure 2. Higher temperatures cause a greater expansion of the plume 
gases through the Ideal Gas Law corresponding to lighter, less dense gases resulting 
in a stronger buoyancy force. The stronger the force of buoyancy the faster plume 
gases will be driven up and away from colder ambient gases. In short, the lower the 
radiation loss the faster the plume gases will rise. Hence, the zero percent radiation 
loss curve is found above all curves depicting a radiation loss, Figures 4 and 5. 

A mistaking assertion of the radiation curves shows several of these curves 
spaced unequally amongst themselves. This is a manifestation of rounding errors and 
not a part of the theory or other underlying assumptions. These curves should be 


spaced equidistantly. 


4.7 Introduction to Near Field Entrainment and Flame Height 

In Chapter 3, the equations for entrainment were modified to develop 
equations for flame height. However, the process is reversed when fitting empirical 
constants to those equations. Entrainment correlations are developed using 
previously determined flame height constants because flame height is easier to 


measure and more accurately quantified than entrainment. 
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In addition, empirical constants for the infinite line are determined with 
knowledge of their axisymmetric counterparts because axisymmetric data is far more 
prevalent and accepted/established than that of infinite line data. More exactly, both 
the profile and entrainment constants from the axisymmetric flame height corrrelation 
are used to determine the correlation for infinite line flame height and entrainment. 
From the theory developed in chapter three, the profile constants for axisymmetric 
entrainment and flame height should equal infinite line constants for entrainment and 


flame height. A flow chart depicting this methodology can be found on page 58. 


4.8 Axisymmetric Flame Height Correlations 

The empirical constants for axisymmetric flame height were developed from 
the data on Figure 6. The region enclosed by the solid black line in the center of the 
graph represents a range of data compiled by Zukoski !!°). Where large fuel source 
diameters have accumulated at the top of the region and small diameters have 
accumulated at the bottom of the enclosed region. Obviously, there is a diameter 
effect exhibited by the data. Although we have included diameter as a term in our 
dimensionless parameters, we are still unaware exactly how and to the extent with 
which it affects our correlations 


Recall the axisymmetric, empirical, flame height correlation from Chapter 3. 
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Figure 6: Axisymmetric Smooth Fit of Dimensionless Flame Height, Z/D, vs. Energy Release 
Rate, Q,* 
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In order to solve for C and C2 two points, not shown on the graph, are taken 
from the enclosed region. The first point is taken from the top of the data range 
corresponding to tall flames. The second point is taken from the bottom of the data 
corresponding to small flames. In an attempt to get the best correlation numerous 


data points were fitted to equation (58). The points below represent the best fit. 


Table 4: Axisymmetric Flame Height Data Fit Points 


a a 
D 


Solving the empirical equation, (58), for C and C2 yields: 


€.='0.152 
Cie O35) 
C, = 0.00590 (from the relation above) 


Substituting these constants into equation 3-(58) yields: 


3/2 az 1/2 7 Z 
QO; = 0.00590 ish sam in al te 
(l-z,)\ D D 


The solid black line on Figure 6 represents this correlation. It can be shown that the 


correlation fits the powers mapped out by Zukoski !?°! shown at the top of the graph. 
For a small flame height to diameter ratio corresponding to a weak energy release rate 
Zukoski finds that flame height over diameter varies with dimensionless energy 
release rate to the second power. Our correlation exhibits the same behavior. If the 
ratio of flame height to diameter is small 1 becomes the dominant term in the second 


set of parentheses. Thus the entire term turns into one and energy release rate varies 
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with flame height to the 4 or flame height varies with energy release rate to the 2"° 
power as was depicted by Zukoski °!, On the other extreme, for a large flame height 
over diameter ratio corresponding to a large energy release rate Zukoski finds that 
flame height varies with energy release rate to the 2/5 power. In our correlation, a 
large Z¢/D will overwhelm the value of 1 found in the second set of parentheses. 
Thus, the simplified equation will yield Z/D°? = Q* or Q*”” = Zp/D as was depicted 
by Zukoski. Between the two extremes, Zukoski finds Z/D = Q**” but says 
realistically the exponent could vary anywhere between 2 and 2/5. A value of 2/3 
should be used with caution. 


[14] 


Also placed on the graph is a correlation by Becker and Liang *’ whose data 


follow the slope of a line exhibiting laminar behavior equal to On . The dashed line 
found in the enclosed region is a well known purely empirical correlation by 
Heskestad °*!. It was placed on the graph as a comparison to our equation. 
Heskestad did not design his correlation to fit the entire range of data that we have 
plotted, but rather used a subset of this range. Heskestad reports that for low Z)/D 
ratios smaller than about 0.5 there is a transition from coherent flaming to distributed 
flamelets, and that data in this area and lower is very scattered. Therefore, 
correlations in this area are subject to a large degree of freedom resulting in a 
substantial difference in the rate of decay between the two equations. It is not 


obvious which correlation is better in this regime. 


4.8.1 Radiation Effects 
Figure 7 shows how our correlation for flame height, presented in figure 6, is affected 
by radiation. Similar to the other graphs portraying radiation we have incorporated a 
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Figure 7: Axisymmetric Smooth Fit Radiation Curves of Z,/D vs. Q,* 
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range of radiation curves from no radiation loss to 60 % radiation loss with 15 % 
intervals between these two extremes. As can be seen from Figure 7, the larger the 
radiation loss the taller the flames. This phenomenon can be explained through an 
extension of the theory presented for the effect of radiation on velocity, Figures 4 and 
5. Recall that less radiation loss results in higher plume temperatures which 
corresponds directly to larger velocities. Flame height is determined indirectly by the 
magnitude of this vertical velocity, w. Remember also from earlier theory that 
entrainment is a function of centerline velocity times a constant. Hence, the larger the 
vertical velocity the faster new air needs to be entrained into the plume to make up for 
this lost velocity, (see the continuity equation). As was shown earlier, flame height is 
proportional to entrainment times the heat of combustion of air, AH,. From here out, 
we will examine the effects of radiation on flame height by looking at the two 
extreme cases; no radiation loss and 60 % radiation loss. 

No radiation loss results in a high temperature plume with a large centerline velocity 
that is capable of entraining large amounts or air in a short distance. Because the 
plume has entrained a large amount of air combustible fuel vapors need not travel 
vary far to find the oxidizer necessary for combustion. Thus, little or no radiation 
loss results in short flames. On the other hand, a plume which has lost a lot of its 
energy to radiation will have low temperatures that produce a weak force of buoyancy 
incapable of driving hot combustible gases upwards at a large velocity. Thus, the 
plume has a weak entrainment rate and combustible gases are required to travel a 


greater distance in search of the oxidizer necessary for their combustion. In 
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conclusion, large radiation losses result in tall flames and small radiation losses result 


in short flames. 


4.9 Axisymmetric Entrainment 

Recall from section 4.7, (Introduction), that although flame height equations 
were developed from entrainment equations the process is reversed when fitting 
empirical constants to entrainment equations. C2 which physically represents the 
expansion of the plume through the entrainment constant, a, should remain the same 
for both of these equations. However, the profile constants C or C; will not be the 
same due to a partial addition of the ‘¥ term in flame height. Thus, both correlations 
will grow the same way with height but will acquire different profile constants as the 
result of different powers of VY. 

Figure 8 shows our correlation for entrainment in the flame zone against data 
compiled by Zukoski ?°!, Recall the axisymmetric, empirical, entrainment equation 


from section 3.9.2: 
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a 8 
where C, = —=;C, =-—a@ 
1 2/3 2 5 


The correlation was determined by keeping the same entrainment constant, C2, 


developed from flame height data and solving for C,./'¥ _ Solving for C, ./¥ is 
easily done by picking:a specific entrainment point from Figure 8; inserting its values 


into the entrainment equation; and then solving for the constant. 
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From Figure 8 the point selected 1s: 
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Recall C2 = 0.357 from the flame height equation. 


Inserting these values into the entrainment equation yields: 


0.536 = C, VP (1.8)"?[1 + 0.357(1.8) 


C,v¥ =0.148 


Through the definition of C, ./* .n was found to equal 9.6. 


AH 
c= c.f “) 
cl, 
a 0.152n 
C,V¥ = No = 0.148 
c,T,sn 1.01(293) 
n=957=9.6 


21] applied an improved version of the k-e-g model of turbulence to 


Tamanini 
buoyancy controlled turbulent diffusion flames issuing from a variety of fuel arrays 
and found that at the flame tip a plume entrains 9 times the oxygen stoichiometric 
requirement. Commonly accepted values of n range from 5 — 15, Delichatsios "°! 
Figure 9 shows entrainment in both the far and near fields on the same graph 
by putting far field entrainment correlations in near field terms. Recall from section 


3.5.1, (Far Field Entrainment for the Axisymmetric Source): 
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To convert far field entrainment into near field terms multiply both sides of 


the equation by Z°” and divide both sides of the equation by D°”. 


3/6 ps 5/3 
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pAgDD? 2 o 


Note- From section 3.5, (Far Field Theory for the Axisymmetric Source), Cy= 


3.87 and C= 0.118 for ¥,-= 20 %. Therefore 2C,C; = 0.169. 


Thus, it is possible to represent entrainment in both the far and near fields on 
one graph. As before, the solid black line represents entrainment in the combustion 
zone, whereas the dashed spurs coming off this line model far field entrainment for 
fire plumes of different flame heights. The graph is read by starting on the solid line 
representing entrainment in the flame zone and branching off onto one of the far field 
spurs at a desired flame height. 

At z/D between 3 and 3.5 both regimes of the fire plume entrain about the 
same amount of air. However, for weak fire plumes exhibiting small flame heights, 
the far field entrains considerably less air than the near field. The opposite 
phenomenon is observed for high energy fires yielding flames larger than z/D = 3.5. 
These large flames result in a smoke plume that 1s capable of entraining larger 
amounts of fluid than its issuing near field. The gaps between various far field 
correlations and the solid near field line can not be attributed to the length of the 
intermittent zone. They are merely a manifestation of the near field correlation being 


incapable of describing entrainment in the far field and vice-versa. 
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4.10 Infinite Line Flame Height 

Figure 10 shows flame height data for the line. The two long lines extending 
from the upper right of the graph to the lower left are linear fits of the data. Notice 
that there is a large gap between these two linear fits ranging from 4.3Z** to 2Z**. 
This discrepancy could be due to individual researcher’s varying definitions of flame 
height itself. More specifically, how they dealt with the intermittent zone defined as 
the fraction of time during which at least part of the flame lies above a horizontal | 
plane located at elevation z above the burner ?°!, For example, did they define the 
end of the flame zone to be where flame exists 10 % of the time, 50 % of the time, or 
90 % of the time? Definition of this zone is extremely important since it can occupy 
anywhere from 40-60% of the total flame height. In most cases, authors did not 
report how they defined flame height. 

Even if two authors used the same definition of flame height it is still quite 
possible to get offset values based on their measurement techniques. Several 
examples of popular measurement techniques include (1) infrared-imaged pictures, 
(2) judgments by eyes, (3) usual photographs, (4) video tape recordings 03] All of 
these techniques could yield different values of flame height assuming the same 
source and flame height definition. For example, eye-averaged results are 
consistently 10 -15 % higher than flame heights determined from video pictures 
assuming 50 % intermittency [20] A detailed explanation of measurement techniques 
and definition of flame height made available by individual authors, where possible, 


can be found in the Appendix A. 
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Figure 10: Dimensionless Flame Height for the Infinite Line 
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As was the case with the data presented by Zukoski in Figure 6 for 
axisymmetric flame height, sources having large widths or diameters tend to 
conglomerate at the top of the data range while those with small diameters or widths 
are found near the bottom of the range. In other words, for a given Z,/D “large” 
sources require less energy release rate per diameter or width, Z**/D, than “small” 
sources. Thus, diameters or widths will increase when travelling along a constant 
Z*/D from right to left. This effect is not that important for large Z*/D ratios found on 
the upper right section of the graph where Z** is the dominant term and D has a 
minimal effect. However, for short flames found on the lower left section of the 
graph the width, D, is the dominant parameter and this effect is extremely important 
when determining whether or not data is laminar or turbulent. 

One of the problems encountered when developing our empirical correlations 
was the possibility of encountering laminar flame height data. Our theoretical 
correlations were developed for purely turbulent data in order to take advantage of all 
the assumptions listed at the beginning of the theory section. If laminar data were to 
contaminate the data pool it could throw off our empirical constants because 
assumptions (1) — (3) would no longer be valid. 

If laminar data is present it will more than likely be found on the right side of 


the data range where the widths are predominantly small. Fluid flow is a function of 


2 


the Reynolds Number, Re = vile . For a given fluid flow and viscosity the smaller 


the width the better the chance of encountering laminar data. Roper '!*! details how 
to predict laminar flames based on nozzle or fuel source size and developed a 


correlation for laminar data which has Z**/D raised to the 2/3 power. Cox "! used 


Al 


this correlation to fit some of his own data which he specifies ie be laminar. Cox’s 
laminar fit, on Figure 10, is the first line all the way to the right of the others that 
comes off of the two linear fits discussed earlier. The other lines are drawn with 
different coefficients but the same slope to pick up other possible laminar data. Thus 
laminar data, if present at all, would be found on the nght side of the data range. For 
this reason our correlation, found on Figure 11, was placed on the upper left of the 
data range. 

Figure 11 shows our correlation of flame height for a line source. The 
correlation is developed by using the same profile constant, C, from axisymmetric 
flame height and solving for the entrainment constant by choosing a “middle of the 
road” data point from Figure 11. Rectangular theory, which will be presented in more 
detail later, says that the profile constant should be the same for both of these extreme 


Cases. 


General rectangular flame height equation: 


S J w!2 (Z ae? Z A\Z 
a Oe eh PCA 1+Cu( Sf 
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The two subscripts, a and b, for the C2 constant are required because the rate 
of entrainment on two different size lengths is not necessarily the same. Where a 
represents the rate of entrainment on the short side and b the long side. In addition, A 


represents the short side and B represents the long side. 


For the axisymmetric case A=B: 
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Figure 11: Dimensionless Flame Height Equation for the Infinite Line 
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As expected, this is the same equation developed for axisymmetric flame 
height. Both C3 and Cz should match axisymmetric “fitted” values. Note that if the 
two sides of the rectangle are equivalent, as in the axisymmetric case, C2, will be the 
same constant as C2». As a result, the multiplication of these two terms will simply be 
denoted C2. In addition, because there is no longer a difference in length between A 


and B the two terms have both been denoted D. 


For the infinite line source A<<B: 


1/2 
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The third set of parentheses on the RHS of the equation reduces to one. In 
addition, the A*” term found on the LHS of the rectangular flame height equation 


becomes minimal when compared to B and thus drops out. Therefore, Q*_, turns into 
O,. 


Comparing the two equations, (A=B and A<<B), it is obvious that the profile 
constant, C3, should be the same for both the axisymmetric and infinite line flame 
height correlations. Recall that we have used the axisymmetric constant to develop 
our linear flame height correlation rather than vice-versa because we felt more 
comfortable with the data available for this source. Note that the entrainment 


constants for the two equations are not the same. 
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To solve for the entrainment constant C2 we set C3 equal to 0.152 from the 
axisymmetric flame height correlation and pick a point from Figure 11. The chosen 
point is located in the center of the top cluster of data. Reasons for picking the point 
here were covered earlier. Recall that the chances of correlating only turbulent data 


are the best at the top of the data range. 


Table 5: Infinite Line Flame Height Data Fit Point 
en 
D 


Solving for C2 knowing C3 = 0.152: 


4.75 = 0.152(10)’*(1+C,(10)) 
9.88 =1+10C, 
C, = 0.888 


Plugging C2 and the axisymmetric, flame height, profile constant into the infinite line, 


flame height equation yields: 


= 3/2 Z 1/2 Z 
QO = 0.00590 fe | fs \Bh 140s 7 
(l-z,)\ D D 


4.11 Infinite Line Entrainment 


Figure 12 shows two possible near field entrainment equations for the infinite 
line. The bottom correlation was developed by using the same profile constant 
developed for the axisymmetric entrainment equation and the same entrainment 
constant from the linear flame height equation, Figure 11. The axisymmetric, 
entrainment, profile constant was put into this equation for the same reasons outlined 
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above for substituting the axisymmetric, flame height, profile constant into the linear 
flame height correlation. 

The top correlation was developed by keeping the same entrainment 
coefficient but determining a new profile constant based on one data point extracted 
from Cox ”!, Unfortunately, we were forced to more or less fabricate this data point 
because to the best of our knowledge no other reliable data was available. The first 
step in creating our mock entrainment point was to take the lowest point on a graph of 
mass flux in a thermal plume by Cox !. Our hopes are that this point is at the 


flame’s tip or just slightly above the flame. 


™ — 0.0055 


I 


v4 
1 


To get a ballpark value for Q; because none were explicitly reported, 100 kW 
was chosen which corresponds to the average visible flame height for turbulent 


flames on Figure 2 of Cox’s paper. 


Thus: 
m, = 0.0055(100) = 0.5548/ 


z = 0.03(100)”"* = 0.647m 


Making the terms dimensionless yields: 


m’ 0.55 
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Inserting the data point into the linear entrainment equation along with the 


value of C2 developed from our flame height correlation and solving for C;: 


aso 5) (He(5) 


78.9 = C,(2.615)V43.1(1 + 0.888(43.1)) 
C, =0.117 


Hopefully, these two correlations found on two totally different premises 
represent the upper and lower bounds of the uncertainty for our entrainment 
correlation. No data exist in the flame zone to verify either of the correlations. 

The process of converting far field infinite line entrainment into near field 
terms is detailed below. Recall the equation for far field entrainment, (16), from 


section 3.4.1: 
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Convert far field entrainment into near field terms by multiplying both sides of the 


equation by Z°” and dividing both sides of the equation by Dag 
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Applying these correlations to the near field, as was done in the axisymmetric 
case, makes little sense because we are not sure whether the near field correlations are 


reasonable. 
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Chapter 5 Empirical Constants 


The five tables below are a compilation of all of the theoretical and empirical 


constants developed thus far. Several key points will be noted. 


Table 6: Far Field Correlations 


Dimensionless Axisymmetric Infinite Line 
Variable 
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Theoretically, if the entrainment coefficient, a, is the same for both the 


combusting zone and the smoke plume, Tables 6 and 5 respectively, the plume will 
expand at a faster rate in the far field for both sources. Yokoi °! documents this 
behavior in Figure 13. (The empirical solutions do not support the constancy of « as 
it is much higher in the combusting region.) The assumption of a constant « and 


similarity profiles close to the fuel source can not be correct. 


80 


Table 7: Near Field Correlations 


Dimensionless Axisymmetric 
Variable 
4 4a 
—a —> 
E | er k 


Figure 3.2 
Horizontal distributions of temperature above circular 
heat source. 


Xre=99 Continuous heat source 
Or. =20 Discontinuous heat source 


Figure 3.3 
Diagrammatical temperature distri- 
butions near a circular heat source. 
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Table 8: Far Field Entrainment Correlations 


Infinite Line 
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Table 9: Flame Height Correlations 


Axisymmetric Infinite Line 
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Further justification for using the same empirical profile constant developed from 


axisymmetric data for the infinite line stems from the fact that the theoretical 


constants, aa and oa are almost equal in the above equations for O7, and Q>. 
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Table 10: Near Field Entrainment Correlations 


Infinite Line 
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Note — The two empirical correlations in the infinite line column represent the upper 


and lower limits of near field entrainment. 
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Chapter 6 Rectangular Theory 
6.1 Introduction 

Having developed empirical correlations for the infinite line and axisymmetric 
source it should be possible to determine correlations for source geometries between 
these two extremes. More specifically, this section utilizes the knowledge gained in 
previous chapters to develop correlations for temperature, entrainment and fants 
height for rectangular sources. The following studies were not exhaustive nor meant 
to be thorough. Their purpose for being was to illustrate the possibility of subtending 
simplified 2-D theory to a 3-D problem. Hopefully, the following discussion will 


spark interest in this possibility. 


6.2 Assumptions 
All of the assumptions listed in section 3.2 are carried over into this chapter. In 


addition, we have added two more. 


(1) Plume cross section maintains rectangular shape. 
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The implication of the first assumption is that the plume will keep its original 
geometry by maintaining crisp comers as it ascends. In reality, the plume will 
become more and more rounded as it rises away from its source. The hazards 
appreciated by incorporating assumption (1) are minor. The second assumption was 
actually proven for both the axisymmetric and infinite line source, and therefore 


assumed to be true here. 


In addition, we will utilize the following Gaussian profiles and set B=1: 
o w yy (=) 
— = ——_ = _-|— ex -_-|j— (2) 
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6.3 Far Field Non-Combusting Plume Temperature 


Assume an appropriate point source energy equation: 


O(- 7,)=4]| p.c,w(T - 7, Jayde (3) 


Integrating this energy equation over % of the source area and multiplying by 4 yields 


the amount of energy liberated over the entire plume cross section. 


Substituting the above Gaussian profiles for vertical velocity and temperature into 


equation (3) yields: 
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Using the following relationship, i} e* adk= “e , the integrals can be solved: 
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Recall equation (1) from section 6.2: 
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Substituting equation (1) into equation (4) yields: 


O(l- z,) ="C.® 3/2 51/2 ob 


Plplov 2 a 
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Note, Q", becomes Q* or Q7 in either the axisymmetric or infinite line limit 


respectively. 
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The above equation is graphed on Figure 14 as y = Cx”. 
Where: 


y=, 


ee (stad) (8) 
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Data points of various aspect ratios ranging from 0.2 to 1.0 were taken from 
Hasemi '"!. In order to convert Hasemi’s data points into our format a value for the 


entrainment constant, C2, had to be assumed. Far away from the source, 


[3] 
F = | , the rectangular plume will resort to point source axisymmetric 


behavior. Therefore, C2 was assumed to equal oa which was taken from 


axisymmetric plume theory developed earlier in section 3.5. The entrainment 
coefficient, @ , was taken as 0.11 from Zukoski Sod However, our value of 0.098 
could just have easily been used. Recall that the difference between these two values 
results from the inclusion of a 20 % radiation loss when determining our value. 

In the combustion region, temperature remains constant. The value on Figure 
13 falls below temperatures for the infinite line and above temperatures for the 
axisymmetric source. 

In the far field, a best fit of the representative data yields a coefficient of 2.3 


which is slightly higher than the coefficient calculated by manipulation of equation 
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Figure 14: Dimensionless Temperature vs. Height for a Rectangular Source 
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(7). As B approaches A in the far field equation (7) reduces to an axisymmetric 
expression for far field temperature. This is rather fortuitous because we can then use 
the far field axisymmetric constant for temperature, Cr, to determine C, found in the 


profile constant of equation (7). The process is illustrated below. 


Setting A = B into equation (7) yields: 


5 bx.) OLa) 


“ (9) 
(nc, KG yy" (2C, iy 
Setting A = B into equation (6) gives: 
AEE: 
ae (=| (10) 


Substituting equation (10) into equation (9) gives: 


ans fey 


C; can now be solved by setting Cr equal to 9.115 from Yokoi "!; C, equal to 6/5a 
from axisymmetric far field theory; a equal to 0.11 from Zukoski ”°!; and y, equal to 
30 %. The constant, C, in equation (8) can now be solved by substituting the value 
for C;. C is found to equal 1.54; a value slightly lower than the fitted value of 2.31. 
The discrepancy between the calculated constant and the fitted constant could be 


attributed to the fact that B was assumed equal to one here in the far field. Because a 


= 0.11 was determined from B + I provisions should have been made to include f in 
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our calculations, or a new @ should have been calculated for B = 1. If either of these 


methods were carried out our calculated C should fall a little closer to the fitted value. 


6.4 Far Field Entrainment 


An expression for rectangular entrainment can be assumed: 


m, = 4l| p,wdydx 


Substituting Gaussian profiles into the above equation yields: 


cree Dae 


m, = 4p,w,,ab || , 


Solving the integrals gives: 


vn Vr 


m,=4p,w,,a Leas 


mM, =1p,W,,ab 


6.5 Flame Height 


(11) 


(12) 


As before, an equation for flame height can be developed by multiplying a 


near field entrainment equation, which includes provisions for source geometry, by a 


parameter containing AH,. 


First, the general entrainment equation, (12), is modified through equation (1): 


m,=7p,C, lg ae I zab 
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(13) 


From equation (5) for a and b and the solution for dimensionless temperature 
in the combustion zone developed below, an expression for near field entrainment can 
be ascertained. 


Assume an appropriate energy equation in the combustion zone: 


A ST (-7,) 


pac, = [Jw(t - T, Ybedy = (1 - x, )ff O"dedy = (4 (14) 


Note that the relation between energy release rate and air entrainment was shown 
earlier in section 3.7.1 equation (33) and section 3.7.2 equation (42). 

In order to solve for ® integrate the LHS with the appropriate Gaussian 
profiles and ignore effects at the source, z= 0. Integrating the entrainment derivative 


on the RHS and substituting from equation (12) yields: 


“ie \t-2,) 


c.w.® T ——=—p 2w_ab 
P©, 4 af m 


0-7 atl Ji-z)=2¥ (15) 


Substituting the values for a and b from equation (5) and ® from equation (15) into 


equation (13) yields: 


m,=7p,C, Jez aan 46c2\F +C,2] 


Simplifying yields an equation for near field entrainment: 
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sachet tad iba ic (=) (i+2c, 2)i+2¢, 2) 
PaizAck B OE ey, A B 


c 1/2 
m v4 
reas = c= c Ge at 28 <) (16) 
| 
4 


Where C, = 


Since m,(Z,)=nsm, = 


, the near field entrainment equation is easily 


converted into flame height. 


ie, 
~ oa) (r20, “tfi+20,4 =} ae 
pice Re Pa nsc 1, 


Orn co (% 24) [is 20 “tl r+2c, (4) (17) 


Because the profile constant, C3, was set at 0.00590 for both the infinite line 


and axisymmetric sources we will use the same value for our rectangular solutions. 
The entrainment constant, C2, was found by fitting the equation to a point in the 


center of the data range. 


Table 11: Rectangular Flame Height Data Fit Point 


Geometry Oe Z/A 


(0.2 x 1.0m 
Note that this point’s geometry should have no “lasting” effect on the flame height 


correlation thereby favoring sources of (0.2 x 1.0 m). 
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The data used to fit this correlation was taken from Hasemi ") and includes 
aspect ratios ranging from 0.1 to 1. The data was plotted directly onto Figure 15, thus 


avoiding any assumptions or manipulation. 
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Figure 15: Dimensionless Flame Height for Rectangular Sources 
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Chapter 7 Conclusion 

Empirical correlations were developed for the infinite line for temperature, 
velocity, and plume width in both the far and near fields through dimensionless plots. 
Correlations were also developed for the axisymmetric source for all of the above 
characteristics. However, due to the acceptance of previous correlations and data 
comparisons of these correlations were not made on dimensionless plots. 

Point and line source theory was modified with a diameter correction to 
determine flame height and near field entrainment correlations for the axisymmetric 
and infinite line source respectively. For the axisymmetric source, results were 
consistent with data and n was found to equal 9.6. For the line source an empirical 
correlation for flame height was developed. However, no entrainment data was 
available for this source making it impossible to develop an exact entrainment 
correlation. As a result, two entrainment correlations were created that hopefully 
embody the upper and lower bounds of a “true” entrainment equation. 

The theory and empirical correlations developed for axisymmetric and infinite 
line solutions were applied to rectangular sources in a fairly rough manner. Future 
research should center around refining this theory and developing more data for 
rectangular sources. In addition, entrainment data for the infinite line source would 


prove extremely valuable. 
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Appendix A Description of Experimental 
Apparatus and Procedures 
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Natural Convection Above a Line Fire 
(Shao-Lin Lee) 


e Temperature 
e Flame Height 


Expenmental Apparatus 


The burner is 78” long by .564” wide and .282” deep. A 1/4” diameter copper circular 
tube is soldered onto both sides of the length of the bumer. The tube is filled with 
water and serves as a water cooling jacket. Fuel is fed into the bottom of the burner 
through six equally separated ports along the length of the burner. Such a 
configuration in conjunction with the “Steady Fuel Supply System” allows for equal 
filling and combustion of fuel in the channel burner. Both methyl alcohol and acetone 
are used as fuels, where the first fuel is used to produce a flame with little heat loss 


due to radiation and the latter produces the opposite effect. 


The burner is housed in a “combustion cage” made out of wire mesh with dimensions 
of 4 feet wide, 8 feet long, 8 feet high and 3 feet above the laboratory floor. The base 
of the combustion cage is made out of 1/4” thick asbestos board and rests on a table. 
On either side of the channel burner, which runs down the length of the combustion 
cage, lie a 1/4” thick 18” wide piece of asbestos board. The asbestos board lies flat on 
the base of the combustion cage and is the exact length of the channel burner. Against 
each end of the side plates stands a vertical asbestos end plate 4 feet wide and 8 feet 
tall which serves to cut off the entrainment of air from the ends. At times the end 


plates are wrapped with aluminum foil to alter the effects of radiation. 
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Temperature measurements are made with a 40 inch piece of pure nickel resistance 
wire 0.002 inches in diameter. The resistance wire is suspended parallel to the channel 
burner between two metal poles and is stretched slightly to account for thermal 


elongation during the experiment. 


Expenmental Procedures 


The combustion cage was centered in the middle of a large room measuring 
approximately 40 feet long, 25 feet wide, and 30 feet high. To insure a minimal 
amount of disturbance in the ambient “8 all detectable leaks around doors and 
windows for example were sealed with adhesive tape. In addition, local heating of 
certain areas by radiation from the sun was reduced by employing window curtains. 
Before each experiment at least one half hour was allowed for the air in the room to 


equilibrate. 
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Flame Height from Rectangular Fire Sources 
(Sugawa, Satoh, & Oka) 
e Flame Height 


Experimental Apparatus 

The gas diffusion burner contained fine sand as a diffuser in a rectangular stainless 
steel vessel. Propane with a monitored flow rate was used in the vessel. The authors 
report aspect ratios of 1:6, 1:8.8, 1:12, 1:17.8, 1:40, 1:60 however, no exact burner 
dimensions were given. In later tests, burner dimensions included 1 cm X 40 cm, 1 cm 
X 60 cm, 2 cm X 60 cm and 2 cm X 120 cm. Most likely the burners used in the first 
set of tests were of the same order of magnitude as those used in the second and third 
set of tests. The authors concluded that flame height is a function of HRR and 


practically independent of fuel and aspect ratio. 


Experimental Procedures 


Three minutes were allowed for flame heights to reach a quasi-steady state at which 
the flame height was recorded by four different methods. 


Infrared-imaged pictures 
judgments by eyes 
photographs 

video tape recordings 


However, the data from the video tapes was primarily used to report flame heights in 


order to eliminate bias inherent in the other three methods. 


Measurements of flame height issued form a single line burner were taken with and 
without both walls and floors. The data reported on my flame length graph is only 


from a burner on the floor and out in the open. 
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Fire Induced Flow in a Vertical Clean Room... 
(Sugawa, Oka, Hotta) 
e Temperature 
e Velocity 


Experimental Apparatus 

All of the experiments were performed in a “clean room” with dimensions of 9m long, 
6.4m wide, and 3.2m high. A 22 cm/sec downward vertical flow was achieved as a 
direct result of the ventilation system which changed the air in the room 200 times/hr. 
The line burner was positioned just off-center at a height 78 cm above the free access 


floor. 


The burner itself was 1cm wide by 85 cm long. Methanol was supplied to the burner 
at a constant rate of .37-.4 ml/sec. The delivery procedure was not documented nor is 


it deemed to be important. 


Experimental Procedures 


Combustion of the fuel was kept constant for about one hour before any data was 
recorded. Thirty three 0.32 mm diameter alumel-chromel thermocouples were spaced 
equally at 10 cm across a movable wire approximately three meters in length. The 
wire was moved horizontally and vertically so as to measures temperatures at different 
points of interest. Thermocouples were maintained at each position for approximately 
two to three minutes. Velocities were measured by a bi-directional tube measuring 
dynamic pressure which was also moved both horizontally and vertically. A smoke 


wire system was used to visualize where velocities were almost zero and indicate a 
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balance between upward and downward flows. The wire was also used to indicate the 


size and position of eddies and flow patterns. 


“Heat flux was calculated based on the consumption rate of the fuel and its heat of 
combustion at 25 °C to gaseous CO, and H2O. The estimated heat flux was about 6 
kW/m’. This is the only HRR given by Sugawa, therefore I followed the author’s 
advice using an average value of the flow rates given and a HRR from Tewarson’s 


chapter in the SFPE Handbook. 


Velocity measurements are taken from approximately 0.2 -1.8 m above the free access 
floor. When extracting these velocity values careful attention must be used because 


the graph starts its datum at the free access floor not at the fire source level. 
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Fuel Shape Effect on the Deterministic Properties of Turbulent Diffusion Flames 
(Hasemi, Nishihata) 
e Temperature 
e Flame Height 


Experimental Apparatus 

The experiments were conducted using a 0.1 by 1.0 m and a 0.2 by 1.0 m rectangular 
porous burner with propane as the fuel. The 0.2 by 1.0 m burner could be sealed off 
into five equally sized square burners thereby making aspect ratios of 1,2,3,4,5, and 10 
available. The burners were filled with ceramic beads of approximately 5 mm in 
diameter. The heat input ranged from 20 kW to 300 kW corresponding to a Q* moa 


between 0.3 and 5. 


The burner was placed 0.50 m above the ground in a large room measuring 27m by 20 


m in plan and 27m high. 


Experimental Procedures 


Flame height was monitored by a video camera during each experiment. The reported 
values of flame height are the average of the height of flame tips observed for more 
than three minutes at the intervals of one sec. The flame heights measured during the 


0.1 by 1 m burner experiments range from 0.1 to 1.5 m. 


Temperature measurements were made using 0.20 mm diameter chromel-alumel 
thermocouples. (The time constant of the temperature measurement system is 
approximately ten sec..) The authors report an under estimation of gas temperatures 


ranging from 2-20 % over 300-1,000 °C. These errors are a direct result of radiation 


105 


and conduction losses at thermocouple junctions and soot buildup on the 
thermocouples themselves. To eliminate the effect of accidental sway of the flame on 
temperature measurements, temperatures were measured not only above the burner 
center but at two different points 5 cm apart from the center in the direction of the 
burner’s shorter side. Reported temperature values are the average of the maximum 
temperature at each height at the interval of 10 s over three minutes during which the 
temperature above the burner center was higher than the other two. The temperatures 


measured during the 0.1 by 1 m burner experiments range from 10 to 1000 °C. 
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Report of the Building Research Institute 
Study on Prevention of Fire Spread Caused by Hot Upward Current 
(Yokoi) 
e Temperature 
e Velocity 


Experimental Apparatus 

The burner was made out of tin and measured 100 cm long by 1 cm wide. Alcohol 
was supplied to the burner via one communicating tube. Yokoi never mentions 
exactly what type of alcohol was used for his experiments however Sugawa reports 
that Yokoi used methyl alcohol on pg. 365 of “Fire Induced Flow in a Vertical Clean 


Room...” 


Experimental Procedures 


Yokoi found that slight draughts markedly disturb upward currents. Therefore, he 
conducted his experiments in a dark room used for photo processing. The room 
measured approximately 4.5m x 4.5m x 4.5m. In addition, experiments were only 
conducted on calm days where there was not much wind and only lasted for a 
maximum of two hours, so as to avoid draughts due to the difference in temperature 


between inside and outside the room. 


A temperature compensating type hot wire anemometer measured upward velocity. 
The meter is not suitable for measuring the velocity of gas at temperatures above 100 


16. 


Gas temperatures were measured using a copper-constantan thermocouple having a 


0.2 mm diameter. The thermocouples were connected to a simple slow rotating 
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oscillograph where the average values were recorded by eye for a time period of ten 
minutes. The automatic recording instrument could record no more than six points at 
a time, therefore one of the points was selected as the standard point and the other five 
points, distributed at various points on the same vertical line, were moved in relation 
to the fire source. As temperatures and velocities gradually changed with time 
adjustments were made in analyzing the results of experiments according to the value 


of the temperature at the base point. 


Both temperatures and velocities were measured between the range of 50-140 cm. A 


HRR of approximately 6 kW was given off by the burner. 
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Experimental Investigation of Laminar Free-Convection Flow... 
(Brodowicz, Kierkus) 
e Temperature 
e Velocity 


Experimental Apparatus 
A wire 0.075 mm O.D. and 250 mm long was stretched horizontally and heated by a 
direct electric current passing through it. Heat supply to the wire was measured at 


9.75 kW by means of a wattmeter. 


Experimental Procedures 


The velocity distribution was determined by introducing small dust particles into the 
air around the wire. These particles were then illuminated by a thin light beam and 
their trajectories were recorded on photograph. Velocity components could then be 
determined from the distance between two positions of a dust particle on a 


photograph. 


The temperature distribution was determined by means of the Mach-Zehnder 


interferometer. 


Measurements of velocity and temperature were taken at points 1, 2, 4, and 8 cm 


above the wire. 
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Gravitational Convection from a Boundary Source 
Hunter Rouse, C. S. Yih, and H. W. Humphreys 


e Temperature 
e Velocity 


Experimental Apparatus 

The room was approximately circular with a diameter of 25 feet and a height of 11 
feet. The burner was housed between two parallel walls measuring 4 feet high by 8 
feet long and 4 feet apart. The burner was placed in the middle of a low platform 
extending the length of the walls. The authors never reported a detailed description of 
the burner including dimensions, fuel used, and other relevant intricacies. The only 


available description describes a “...recessed gas burner yielding low, blue flames...” 


Experimental Procedures 


Temperature measurements were made by a copper constantan thermocouple in 
combination with a potentiometer reading to 0.002 millivolt. Velocity measurements 
were obtained from a vane anemometer | 1/4 inch in diameter. Velocities as low as 
0.2 feet per second were measured. Both instruments were mounted alternatively on 


remotely controlled traversing mechanisms. 
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An Experimental Study of Some Line Fires 
Li-Ming Yuan and G. Cox 


e Temperature 


e Velocity 
e Flame Height 


Experimental Apparatus 

Natural gas, 94% Methane, was burned in three different burners including two porous 
refractory burners measuring 15mm x 0.20 m, 15 mm x 0.50 m, and a sandbox burner 
50mm x 0.5m. The fuel supply could be adjusted to achieve theoretical heat release 
rates ranging from 2 to 110 kW. Similar to the setup by Lee and Emmons, 
experiments were conducted in a large laboratory with a fine mesh screen surrounding 
the burner to minimize the effects of ambient air movements. The burner surface was 


elevated 0.7 m above the floor. 


Experimental Procedures 


Flame height measurements were obtained from 3 min samples of video records. 
Temperature measurements were made using 150 micrometer chromel alumel 
thermocouples, and velocity measurements were recorded through 10 mm diametre bi- 
directional pressure probes. Both sets of measurements were taken up to 1.2 m above 
the burner and midway across its width. Similar to flame height measurements, both 
centerline time-mean temperatures and gas velocities where determined over three min 
samples. More specifically, temperature measurements were determined by averaging 


the maximum temperature exhibited by three thermocouples at the same elevation. 


11] 


One thermocouple was placed at the center of the burner with two the other two offset 


50 mm across the burner’s width. 


Mass flux was also determined for the 0.5 m long burners through a hood and duct 


technique. 


Linear Flame Heights for Various Fuels 
F. R. Steward 


e Flame Height 
Steward did not describe any of the experimental procedures or the apparatus used to 


develop the data he used in his paper. 


The Size of Flames from Natural Fires 
P. H. Thomas 


e Flame Height 


Thomas also has limited information regarding experimental procedures and apparatus. 
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Appendix B Data in Chart Format 
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Experimental Investigation of Laminar Free-Convection Flow... 
(Brodowicz, Kierkus) 
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An Experimental Study of Some Line Fires 
Li-Ming Yuan and G. Cox 
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burner dimensions 
© 0.5m tong x 0.015m wide 
6 0.5m tong x 0.050m wide 
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Fig. 3. Reduced centreline velocity with height. 
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Visible flame height L (m) 
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Fig. 2. Visible flame height as a function of Q,. 
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Experimental study of line fires 
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Fig. 4. Centreline temperature rise with reduced height. 
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Fuel Shape Effect on the Deterministic Properties of Turbulent Diffusion Flames 
(Hasemi, Nishihata) 


Flame Height 
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Hasemi (200 kW) | Hasemi (331 kW) 
_AT/T, 
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La SR ee 


Hasemi (6.4 kW) Haseml (34.2 kW)" Haseml (170 kW)* 
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Figure 6 

Maximum excess Remperacugs /xs- 
normalized height(z/Q A) 
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Natural Convection Above a Line Fire 
(Shao-Lin Lee) 


Flame Height (Asbestos/Methanol) 


Se ae 


or | 2 | zp | 
05 78 
44 


75 3.102837 0.00156  2.9796| 2.070643| 0.019133| 1.337954 
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Figure 5-11 Measuroment;s of Average Height of Acetone Plame 
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Figure 5-lt Hearuremonts of Average Ileight of Methyl Alcohol Flames 


131 


_ S15 


_ ACETONE 


WITH ALUMINUM 
FoI. COVERED 
SIDE PLATES 


a 


218 
Ay 


e 

g 

a 

Q 10 

Ww 

-) 

= 

g 

& 

at 

ry) 

Y) 

= Ss © $= 0.0187 gu/Sec-—m 
a 0.0165 7) 
ce) ©0127 °° 
v 0.0/03 
0 0.0070 “ 
+ 0.004] oo 
a 0.00258. ry 
x 0.00157 « 
eo 0-001 34 
® 0030/0 


fe] 10 20 30 +0 50 60 70 
Vertical Distance X, 'n, 


Firers 5-7 Suovancy Weasuremomss for the Flumo above a Flame of Asetons 
with Alumimm Foil Covered Side Plates = 


132 


0 Faoceise gn/me-em) wiry BARE ASBESTOS 
pnd eee " SIDE PLATES 

. © 0. OOOTIE o 

fy & 000158 » WITH ALUMINUM FOIL 
Vv 0. 00074 COVERED SIDE PLATES 

BSS 

w 40 

ss 

9 

Ss 

i @) 

4 30 

3 

~ 

co 

1S) 

© 20 

‘a 

yy) 

cc 

to METHYL ALCOHOL 


° 10 20 30 40 50 7) 70 
; Vertical Distance =X , én. 


Figure 5-3 Buoyancy leasuroments for the Plume above a Vjame of Methyl Alcohol 
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Gravitational Convection from a Boundary Source 
Hunter Rouse, C. S. Yih, and H. W. Humphreys 


| Q@ | Zz [| zz | 22 | dT, | 
| 0.752981|0.01_| 51.80 [0.019303 | 0.05 _| 
| 1.003975] 0.01 | 49.34 0.020269 | 0.05 _ 
0.752981! 0.01 | _75.71 [0.013208 |_0.03_| 
| 1.003975] 0.01 | 65.78 |0.015201 | 0.04 
| 1.039202{ 0.01 | 70.72 _|0.014141 | 0.04 _| 
1.166901] 0.01 | 101.16 0.009885 | 0.02 _| 
| 2.298575] 0.02 | 66.25 __|0.015093 | 0.04 _| 


0.39624 
0.4572 
0.57912 
0.6096 
0.67056 


1.0668 
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HUNTER ROUSE, C. S. YIH, AND H. W. HUMPHREYS 
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Fig. 2. Distribution functions for two-dimensional convection. 
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Fire Induced Flow in a Clean Room with Downward Vertical Laminar Flow 
(Sugawa, Oka, Hotta) 
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Figure 4 Variation of upward 
velocity along the ‘center line. 
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- Yokoi (Temp.) 
: Yokoi (Vel.) 
: Temp. 
: Vel. 
:Vel. without downward flow 


our data 


Temperature(°C), Upward velocity (cm/sec) 


1 5 10 50 100 200 
Height from fire source (cm) 
Figure 9 Variations of temperatures and velocities 


along the center line of the upward flow with and 
without down flow. Yokoi's data are also plotted. 
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Flame Height from Rectangular Fire Sources 
(Sugawa, Satoh, & Oka) 


3.117611 
“7638039 


4.085778 
4.167949 
4.645153 
4.980012 
4.980012 
5.245826 
5.813898 
6.087695 
7.978215 
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y without floor 


g with wall 


Line fire 


0.1 05 1 a 5 10 50 100 

Q rec 
Figure 2 Dependence of flame height in dimensionless form on heat 
release rate in dimensionless form. Data were obtained from a Lise 


fire source using an unconfined rectangular propane gas diffusic= 
burner without wall (0, v) and against wall (s). 
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Linear Flame Heights for Various Fuels 
F. R. Steward 


x-axis 
100 
112.2018 
199.5262 
316.2278 
446.6836 
562.3413 
794.3282 
63.09573 
50.11872 
35.48134 
28.18383 
17.78279 
12.58925 
10 
8.317638 
5.623413 
4.265795| 1621.81 
2.511886| 1359.356 
1.584893 
1.513561 
0.630957 


5. 0.316228] 681.2921} 6.329203 


1.2: 
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0.239883 
0.177828 
0.158489 
0.125893 
0.063096 


Hydrogen) 
x-axis 
0.398107 
0.562341 
1.584893 
1.659587 
2.511886 
3.548134 
3.548134 
4.466836 
7.943282 
11.22018 
12.58925 
19.95262 
62 41.68694 
50.11872 
85. 70.79458 


1.97, 93.32543 
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Figure 3. Flame she to base width ratio versus 
mixing controlled burning rate parameter 
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The Size of Flames from Natural Fires 
P. H. Thomas 
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MODELING PRINCIPLES 


L = 400 (m'_) 7”? c.9.8. UNITS 


——=— X LONG STRIP SOURCE 
(PHOTOGRAPHED ) 
MFT <D<2FT 
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U 10-2 10°? 


a DIMENSIONLESS 
PY) ss 


Fia. 7. Flames from long strips and windows. 
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Report of the Building Research Institute 
Study on Prevention of Fire Spread Caused by Hot Upward Current 


(Yokoi) 


146 


SEE aan a 
Bennie 


Figure 1.7 

Vertical distribution of temperature and upward. 
velocity at the axis of upward current from a 
line heat source. 
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heat source (z cm) 
K 


Height from a fine 


of temperature 


40 |— OMeasured values 


30 
0 20 30 40 60 80 100 
In the case to use temperature co-ordinate (C) 
In the case to use upward vel. co-ordinate (em:sec) 


Table 1.2 Comparison of the calculated values with the measured values of the 
temperature at the central axis of the upward current. (in the case of 
an infininte line heat source) 


' Values of temperature at the 


Height from Measured values Absolute temper- axis of stream, J6m 
a heat of ascending ature of the 

source, 2 velocity, Wm ambient air, 60 calculated j measured 
SO cm 103. cm/sec 290 *K 4G : 44 »°°C 
6c" 105 ® ” 340C«n | 33 0=O 
70 " 105 " 292" 3200” 
8&0 7” 108 ” ” 26 ” | 27 w 
90” 102 ” ” 23 ” : 22 ” 
100 “ | 107.5% wv 20 ” | 21 ” 
lO” | 103 " 18.6" | 18.5 
120 104 " “) 17e\ ee ! 17 ” 
Nelo | 106 ‘J 18.8 146.5 
140 "” | 108 # ” 14.6 * | 16 as 
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1x C Conversions of Data 
Appendix 
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Conversions 
In order to consolidate the data presented by numerous authors onto one single 
graph conversions between the data presented in its onginal form and that used in my 
graphs were often necessary, as would be expected. Listed below are a summary of 
the variables used in the plots; definitions of these variables; and generic numerical 
values used in the variables unless otherwise stated by the individual data sets 
themselves. Note, the notation used in this appendix may differ from the rest of the 


thesis. 


Numerical Values: 

Unless stated other wise ambient values are assigned as follows: 
Specific Heat (c,) 1.01 kJ/kg-K 

Ambient Temp (T,) 293K 

Ambient Density (p.) 1.21 kg/m? 

Gravity (g) 9.81 m/s” 


Variable summary: 


cm (2 =) 


Plot summary: 


Temperature Velocity: Flame Length: 
; L 
Ordinate: AL Ordinate: = Ordinate: —~ 
1, V8z* w 
' Vd 
Abscissa: —_ Abscissa: ses Abscissa: — 
Zz z* Ww 
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Experimental Investigation of Laminar Free-Convection Flow in Air Above Horizontal Wire with 
Constant Heat Flux 
K. Brodowicz and W. T. Kierkus 
Temperature Conversion: 
Data was taken from “Experimental Investigation of Laminar Free-Convection Flow in Air Above 
Horizontal Wire with Constant Heat Flux.” The values were graphed against T - T.. CC) and y (mm). 
Only centerline data was extracted therefore there was no use for the x-axis, y (mm), which represents 
distamce off set from the centerline. Temperature values were taken at arbitrarily chosen levels above 
the wire including 1, 4, and 8 (cm). 
Ordinate: T - T.. CC) 
Absissa: y (mm) 
(1) Divide T - T,. or AT by T, = 19.8 CC) = 293 K 


(2) Determine z*: 


2/3 2/3 
z*= (2 ~ ee) —- 423x107 (m) 
ples 112152 


(3) Convert arbitrarily chosen z (cm) into z (m) and divide by z* (m). 


Velocity Conversion: 


Data was taken from “Experimental Investigation of Laminar Free-Convection Flow in Air 


Above Horizontal Wire with Constant Heat Flux.” The values were graphed against u (cm/sec) and y 
(mm). As before, the x-axis, y (mm), was ignored because only centerline values were used. Velocity 
values were taken at arbitrarily chosen levels including 1, 2, 4, and 8 (cm). In order to graph these 


values against u/ 82 * velocity values were divided by 0.0644 and z was divided by z* = 4.23 x 10° 
4 


(1) Determine z* using the procedure outlined above. 


(2) Convert arbitrarily chosen z (cm) into z (m) and divide by z* (m). 


(3) Determine 4 8z ei 
\gz* = [981 )(423 x10 m) = 0.0644— 
S S 


(4) Convert u (cm/sec) to u (m/sec) and divide by ./2z * or 0.0644 m/sec. 
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An Experimental Study of Some Line Fires 
Lin-Ming Yuan, G. Cox 


Flame Height Conversion: 
Ordinate (Visible Flame Length): L (m) 
Absissa: Q; (kW/m) 


(1) Divide flame length L (m) by the burner width, w (m) 


L, (m) 
0.015 (7m) 


(2) Determine z*: 


: 2/3 Pr 2/3 
zt= [2 = ( g = 0.0334 (mm) 
Pole VE 112152 


(3) Divide z* by the burner width, w (m): 


z* (mm) 
0.015 (m) 


Temerpature Conversion: 
Ordinate: AT (°C) 

Absissa: Z/Q,”* (m-kW) 
(1) Divide AT by T, or 293 K. 
(2) Determine z/z*: 


2/3 


z 2/3 z 
—, x(p.T,c, V8) mi QO? x (112152) 


1 1 
Velocity Conversion: 
Ordinate: w/Q,’” (m-s'-kw"”) 


Absissa: Z/Q,?? (m-kW”) 


(1) Determine ee ; 
gz 
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1/3 


air * (e07 ZO array gir * (12152) a7 ae 


(2) Determine z/z* using the same method outlined above. 
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Fuel Shape Effect on the Deterministic Properties of Turbulent Diffusion Flames 
(Hasemi, Nishihata) 


Flame Height Conversion: 


Data was taken from “Fuel Shape Effect on the Deterministic Properties of Turbulent Diffusion Flames” 
page 279, Figure 3. Like Sugawa’s data, X-axis values had to be raised to the 2/3 power so that they may 
be graphed on the z*/w axis. 


Ordinate: O° 


oq (dimensionless) 


Absissa: LA (dimensionless) 


; where B = burner length and A = burner width. 


ee Soon le 
mod pcT. |g A>"B 


« 5/2 3/2 3/2 
Ma a ee 
mi w?] w2 w 


aa" 


Note that the ordinate is already in the correct form and no conversions are necesary. 


Temperature Conversion: 


Data was taken from “Fuel Shape Effect on the Deterministic Properties of Turbulent Diffusion 
Flames” page 280, Figure 4 and page 282 Figure 6. The values were graphed against © ((C)and 


z/ OF eA. After a little manipulation it was apparent that z/ Or eA equals z/z*. The only 
adjustment required to place these values on my graph was to divide the Y-axis by 293 K. 

Although Hasemi plotted values for 15 different heat release rates only nine of them were 
graphed on the regular temperature graph so as not to over shadow the values listed by the other authors. 
As a general rule the two extreme heat release rates were graphed followed by every other value in 
between so as to get an unbiased distribution. One additional heat release rate was graphed to fill a hole 
not covered by the other values. 


Ordinate: z/O*., *° @ A (dimensionless) 


Absissa: © (°C) 
Os a= pest where B = burner length and A = burner width 
pcT,[gA*"B 
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From above flame height conversion: 
2s Were 


Ons =A 


Zz ! Z 
z* 


ENOCH, = 0A z* 


(1) Divide ©, which equals AT, by T.. 


(2) The x-axis is alreay in the correct format and can be left alone. 
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Natural Convection Above a Line Fire 
Shao-lin Lee 


Flame Height Conversion: 


Data was taken from “Natural Convection Above a Line Fire” pgs. 5-21 and 5-22. 
The flame height was given in inches and was therefore divided by the burner width in 
inches to yield L¢'w. To graph the x-axis values given in terms of a fuel consumption 
rate, (q), on a z*/w axis requires a little more manipulation. First multiply by the 
appropriate AH (chem) for either acetone or methanol, and then multiply by 100 to 
convert cm into m. AH chem) Values were taken from Tewarson’s chapter pg. 3-78 in 
the Second Edition of the SFPE Handbook. Next raiseQ’ to the 2/3 power and 


multiply by .00924 which yields z*. This value is then divided by the width of the 
burner converted into m or .0143 m. 


Ordinate: L (in) 

Abscissa: q (gm/s-in) 

(1) Divide the flame height, L, by the width of the burner 0.564 in to yield L¢w. 
(2) Find QO: 


Q' = qAH, 


AH rote 
; & 


Cn 


AH. = 279. 
& 


Cmmeth 


A quick example best explains this step. 


OU = Ghia (0001ss—£—\( 279) 00em = 43 
= g lm S—m 


(3) Raise 4.3 kJ/s-m to the 2/3 power. 
(4) Multiply the resultant by 0.00924 or- 
Oe . 


92/3 
oa Ss 


(c.c,T.ve) tal 11215222 


(5) Divide z* by the burner width or 0.564 in converted to 0.0143 m to get z*/w. 
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Gravitational Convection from a Boundary Source 
Hunter Rouse, C. S. Yih, and H.W. Humphreys 


Temperature Conversion 


—Ay ie AT IT, ~ AT/T, 
(p(—u/ L)’ eae ea ( Q’ / ZaAZ 


<a 


(1) Convert x(ft) to x(m) 


(2) O/i= C,T,(w/L) look at eq. 12 & convert to SI units 


= 101 oa ao se ) asa #8) a 
kgeK L se ft lb/ 3048m 


= 40472 | 
(3) 
Find 
z+-( Ol ) = Q 4 (2) y" 
Cig 12()10f % ag 9x oa (2 =) eer 


o(Z)-8 


Note that the Y-axis can be left as is because it is already in the correct form. Z/Z*, 
however, was transformed into SI units for ease of operation. 


156 


Temperature conversion: 


Temperature data came from pages 5-6 and 5-15 of Shao-lin Lee’s thesis, “Natural 
Convection Above a Line Fire.” 


Ordinate (Reciprocal of Buoyancy): fs (dimensionless) 
yA 


Ty £0, ape 1s 


1 — oro _—__"o  —__§2__ because... 
Ay s(e,~P) (e,-p) T-f 


Absissa (Vertical Distance): x (in) 
Fuel Flow Rate Variable: q (g/s-cm) 
(1) Convert x (in) to z in m. 


(2) Find z* using the conversions illustrated above for flame height. 


(3) To get =. divide step one by step two. 


(4) = is achieved by taking the inverse of the ordinate. 


oO 
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Fire Induced Flow in a Clean Room with Downward Vertical Flow 
Sugawa 


Temperature Conversion: 


Data was taken from “Fire Induced Flow in a Clean Room with Downward Vertical 
Flow” page 364. The values were graphed against AT (°C) & H (cm). In order to 
graph these values against T/T.. and Z/Z* temperature values were divided by 293 K 
and Z was divided by Z*=.0334. The value for Q used in Z* was found by taking the 
average volumetric methanol fuel flow rate for the entire burner length times the 
density of methanol times the chemical heat of combustion for methanol (AH chem). 
Sugawa never indicated whether the center line temperatures were recorded with or 
without downward flow. However, in my opinion ignorance of this fact is not 
detrimental as the flow was not fast enough to significantly alter temperature readings. 


Ordinate: AT (°C) 
Abscissa: H (cm) 
(1) Divide temperature by T,. 


(2) Calculate Q’: 

Sugawa reported a range of fuel flow rates between .37 and .4 ml. Therefore the 
average of the two flow rates will be used to determine the heat release rate per unit 
length. In addition, the value used for AH, was found in Tewarsons chapter of the 
Second Edition of the SFPE Handbook on pg. ? 


(37.4 


O= Ga x (.793-£] x (isa ¥) = 583 kW 
2 ml gZ 


. 583kW 
'= (eel 


85m 


(3) Determine z*: 
2/3 


5836” 
el Pe | Al 
cas Crear Cary 
1.21—=— |(293 K) 101————_ 
( m* ( ) kg-K vg S 


(4) Convert H (cm) into (m) and divide by z*. 


Velocity conversion: 
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Data was taken from “Fire Induced Flow in a Clean Room with Downward Vertical 
Flow” pages 361 & 364. The values were graphed against u (cm/sec) and H (cm) on 
page 364 and u (m/sec) and H (m) on page 361. In order to graph these values against 
u/(gZ**)"” and Z/Z** velocity values were divided by .572 and Z was divided by 
Z**= 0334. To see how either of these values was attained reference Sugawa 
temperature notes. 

Velocity values were recorded from two different graphs off the same paper under the 
same general conditions. Intuition tells us that these values should be the same, 
however they are slightly different. Unfortunately, I have no understanding of the 
discrepancy. 


Ordinate: u (cm/sec) or (m/sec) 

Abscissa: H (cm) or (m) 

(1) Convert both axis to meter based values if necesary. 

(2) z* has already been computed above under temperature conversions. 


(3) Divide z (m) by z* (m) or .0334 to get z/z*. 


(4) Find fgz*: 


Jgz* = (981 )(0.0334 m) = 0.572 = 
S 


(5) Divide u (m/sec) by /8z * (m/sec) to get : 


faa 
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Flame Height from Rectangular Fire Sources Considering Mixing Factor 
(Sugawa, Satoh, Oka) 


Data was taken from “Flame Height from Rectangular Fire Sources Considering Mixing Factor” page 
438. The values were graphed against L¢D vs. Q*rec. In order to graph these values against z*/w the 
Q*rec values had to be raised to the 2/3 power. 

Ordinate: L/D (dimensionless) 


Absissa: OX. (dimensionless) 


O- cy wi near ti 
ie asa Ties 


; where W = burner length and D = burner width. 


e 5/2 3/2 
%3/2 * 
Ouse eee & 
rec w?2] w? w 
(=) e 2/3 
Oke 
Ww 


vanes 
Ps ee 
Remember: Z -( a >W=1;D=w. 


Note that the ordinate is already in the correct form and no conversions are necesary. 
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Linear Flame Heights for Various Fuels 
F. R. Steward 


Flame Height Converson: 


Data was taken from, “Linear Flame Heights for Various Fuels” on page 177, Figure 3. Neither the 
ordinate nor the absissa were plotted on a logarithmic scale, however, data on both axis was plotted as a 
logarithmic function. The y-axis need only be converted out of its logarithmic function and it is then 
ready to be plotted. The x-axis also needs to be converted however it rquires the removal of many terms 
found in its equation before it is in a suitable form for plotting. 

Values for heat release rate and stoichiometric ratio were found in Archibald Tewarson’s chapter of the 


SFPE Handbook, Second Edition, on pa th exception of the heat release rate for hydrogen. This 
value was not reported by Tewarson and calculated. 
L’\r+op,/p,) @ 
Ordinate: og EM + 04/0.) @ 
W gQ: Pa (1 = @) 
Absissa: log H/W 
Fuel Stoichiometric Heat of combustion Dasdanoe 
mass of air/mass of fuel P= 
req’d for combustion Pair 
Methane 17.2 50,200 0.553 
Propane 15.7 46,000 1.52 
Hydrogen 34.7 120,900 0.0348 
)., was developed through manipulation of the Ideal Gas Law as follows: 
R PM 
P= ar =f *\r— = — 
A i, ae SIRT. 
Prair = M,, ; thus: 
P fuel fuel 
1.0080 44 16 
? =— = 0.0348; p!. ,, =—— Dien =— =0553 
Pow 98.95 Pocstte — 98 95 Porte ~ 9895 
Calcualtion of the heat of combustion for hydrogen: 
1 
H, ete —+ H,O+Q 
kJ 
= H, -H, = 24183—————_ - 0 
eae i mole H,O 
Ne eh. 2 24183 kJ/mole H,O = 24183 kJ ‘ 1 mole H, ~ 1209 4, 
My  1moleH, /mole H,O mole H, 2g g 


Note the best way to illustrate this conversion is through an example. Therefore the methane point 
| 2 
(0/L)\r+@p,/p.) @ 


log 
W* gO: p;(1-@)’ 


= -~(.030 and log H/W = 1.205 will be examined. 
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(1) Determine Lf/D: 
log H/W = 1.205=> H/W = 10'”” = 16.03 


(2) a= acm = co = 0.0921 ;(from manipulation of eq. 11) 
joe 50,200 — 
clr a kg 
(101 a! Jess K\172) 
kg — K. 
er 0.0921 
(r+0/p,'} oO (172+ al 0.0921 


2 -1 
= 1.042 x10” (_ 
kg-—s-—m 


pa(l-@) (981%][ so 200 A) a2 “81 - 0.0921)’ 


Note - 9/, =Pof) 


log (x-axis) = -0.030 => (x — axis) = 10°” = 0.933 


o93a{—_1__) - 8.95x 10° = ae Di 

1.042 x 10° 

eee = (8.95 x 10°)” 

( Q ie 1 


1 2/3 Ps 2/3 
z 8\1/3 
== =(895x10 (__1 __| = 8.947 
W rea T, J/g ea 121x 101 x 293 x 981 


Zz * 
— = 8.947; whereW=w 
w 
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The Size of Flames from Natural Fires 
P.H. Thomas 


Flame Height Conversion 


Data can be found in “The Size of Flames from Natural Fires” on page 850. The experiment was 
referenced from a paper conducted by Thomas, Pickard, and Wraight “for wood fires on effectively 
infinite strips of width D.” In order to graph Thomas’ data on a Z**/D axis his X-axis values were 
multiplied by 43.9 and then the total was raised to the 2/3 power. 43.9 was attained from 
AHc/(To*Cp). 13.0 kj/g was used for the heat of combustion. According to Thomas his first four 
values do not follow the same slope of his latter due to the discrepancy at small values of L/D which 
arise from the difference in geometry between the two situations. 


Ordinate: L/D 


aff 


p.v\gD 


m” m [Fee 


Absisca: (dimensionless) 


wae SS Oe 


p.jgD  p\giD® “\AH.Te, 


Note - Q = mAH, 


baal (5) GE 
| x 
Pelee DV IAH, \\D AH 


Thus: 


2/3 
( ) [ 7 ( ; 
D PovgD cols 


Where a value of 13.0 kj/g was used for AH, because the specie of wood used during the test was not 
reported. 


13 10° 
AH, = ee = 430 
; 
cp, (10197 x)(293K) 
Therefore: 
z* 


a (43.9(x - axis)) a 


The ordinate is already in the correct form and no changes are required. 
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Report of the Building Research Institute 
“Study on the Prevention of Fire Spread Caused by Hot Upward Current” 
Yokoi 


Yokoi never reported a heat release rate, therefore it was assumed that the heat release rate per unit length 
for his burner was equivalent to that used by Sugawa in his paper, ”Fire Induced Flow in a Clean Room 
with Downward Vertical Laminar Flow,” because both burners had the same width and used similar fuels, 
acicohol and Methnol respectively. The value calcualted for Yokoi’s burner was compared to heat release 
rates he reported for axisymmetric burners using an “equal area concept.” The two values are relatively 
close to each other. 


Calculation of heat releae rate using Sugawa’s burner: 


( 37+4 zl ey 
O= S85 (793 £) y C4 = 583 kW 
2 ml g 
583 kW kW 


O' = ——— = 583 — 
lm m 

Note - For more detailed information reference Sugawa’s paper. 
Calculation of heat release rate using equal area concept: 
Area of line burner = 1 cm x 100 cm = 100 cm? 
Interpolating from 
Temperature Conversion: 
Data was taken from “Upward Current from an Infinite Line Fire Source” page 11, Figure 1.7 and page 
25, Table 1.2. The values were graphed against AT (C) & H (cm). In order to graph these values against 
T/T., and Z/Z** temperature values were divided by 290 K and z was divided by z* = .0337. 
Ordinate: AT (°C) 
Absissa: H (cm) 
(1) Convert H (cm) to z (m). 
(2) Divide temperature by 290 K. 


(3) Determine z*: 


: mh 2/3 
zt= — = (£36) = 0.0337 (m) 
Poloc, V8 1110.04 


(4) Divide z (m) by z* (m) or 0.0337 (m). 


Velocity Conversion: 
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Data was taken from “Upward Current from an Infinite Line Fire Source” page 11, Figure 1.7 and page 
25, Table 1.2. The values were graphed against u (cm/sec) and H (cm). In order to graph these values 
against u/ “J gz * and z/z* velocity values were divided by .575 and Z was divided by z* = 0.0337. 
Ordinate: u (cm/sec) 

Absissa: H (cm) 

(1) Convert H (cm) to z (m). 


(2) Divide z (m) by z* (m) or 0.0337 (m) as previously determined in the temperature conversin section. 


(3) Determine Ri gz* 


Jgz* = (oa) (0.0337 m) = 0575—~ 
Ss 


(4) Convert u (cm/sec) to u (m/sec). 


(5) Divide u (m/sec) by gz * (m/sec) or 0.575 (m/sec). 
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